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University  of  Texas  at  Austin. 
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1.  INTRODUCTION 


The  empirical  characterization  of  linear  time-invariant  (LTl)  systems  is  a  subject  that 
is  well-established  in  the  technical  literature.  The  use  of  deterministic  functions  (e.g.,  unit 
impulse  functions,  unit  step  functions,  and  controlled  sinusoids)  for  the  characterization  of 
LTI  systems  is  intimately  related  to  concepts  governing  the  analytical  characterization  of 
these  systems  (e.g.,  the  unit  impulse  response,  the  unit  step  response,  and  the  system 
frequency  response).  The  use  of  non-detenministic  signals  for  the  characterization  of  LTl 
systems  is  intimately  related  to  concepts  governing  the  analytical  characterization  of 
stationary  random  processes  (e.g.,  covariance  functions  and  the  power  spectrum).  In 
either  case,  the  empirical  characterization  of  LTl  systems  proceeds  amid  a  framework  of 
familiar  analytical  concepts. 

The  analytical  characterization  of  linear  time-varying  (LTV)  systems  is  also  well- 
established  in  the  technical  literature.  Appropriate  extensions  of  concepts  governing  the 
analysis  of  time-invariant  systems  provide  useful  insight  into  the  behavior  of  systems  with 
known  time- varying  parameters.  The  empirical  characterization  of  continuous-time  LTV 
systems  has  received  some  attention  in  the  technical  literature,  but  the  subtleties  of  applying 
these  techniques  to  discrete-time  LTV  systems  and  their  practical  limitations  have  not  been 
discussed. 

The  objective  of  the  following  work  is  the  application  of  techniques  for  empirically 
characterizing  discrete-lime  linear  periodically  time- varying  (LPTV)  systems.  Specifically, 
we  wish  to  characterize  certain  LPTV  structures  found  in  multi-rate  quadrature  sampling 
and  quadrature  demodulation  systems  (see  Baugh,  1988).  Chapter  2  presents  methods  for 
the  analysis  of  LTV  systems  in  both  the  time  domain  and  in  the  frequency  domain.  These 
analytical  methods  rely  on  the  concepts  of  the  time-varying  system  response  function  and 
the  system  transmission  function  (see  Zadeh,  1961,  Zadeh  and  Desoer,  1963,  D'Angelo, 
1970,  Claasen  and  Mecklenbrauker,  1982,  and  Crochiere  and  Rabiner,  1983).  Chapter  3 
presents  a  discussion  of  the  analysis  of  random  proces.ses  generated  by  LTV  systems.  The 
evolutionary  spectral  representation  of  a  nonstationary  signal  and  its  relation  to  the 
underlying  LTV  system  model  is  presented.  Techniques  for  the  empirical  characterization 
of  discrete-time  LPTV  systems  are  presented  in  Chapter  4.  These  techniques  are  presented 
within  the  framework  of  concepts  for  the  analytical  characterization  of  LPTV  and  the 
signals  produced  by  such  systems  di.scus.sed  previously.  The  limitations  of  application  of 
these  empirical  techniques  in  practice  are  also  presented.  The  results  of  applying  these 
techniques  to  LRl'V  systems  of  interest  are  pre.scnted  in  Chapter  5. 


2.  REVIEW  OF  THE  THEORY  OF  DISCRETE-TIME  LINEAR  TIME- 
VARYING  SYSTEMS 

The  theoiy  and  analysis  of  discrete  and  continuous  time  linear  time-varying  (LTV) 
systems  has  been  the  subject  of  research  for  some  time  (see  Zadeh,  1961,  Zadeh  and 
Desoer,  1963,  D'Angelo,  1970,  Claasen  and  Mecklenbrauker,  1982,  and  Crochiere  and 
Rabiner,  1983).  Characterizations  of  such  systems  in  both  the  time  domain  using  the 
concept  of  the  time-varying  system  response  function  and  in  the  frequency  domain  using 
various  concepts  related  to  the  system  response  function  have  been  proposed.  In  this 
chapter,  we  will  review  some  of  these  concepts  for  the  characterization  of  discrete-time 
LTV  systems.  For  a  more  thorough  treatment  of  the  subject,  see  D'Angelo,  1970, 
Crochiere  and  Rabiner,  1983,  and  Zadeh  and  Desoer,  1963. 

2.1  TIME-DOMAIN  REPRESENTATIONS  OF  LINEAR  TIME- VARYING 
SYSTEMS 

A  discrete-time  linear  time-varying  system  may  be  described  mathematically  as  a 
linear  mapping  of  input  signals,  x(t),  to  output  signals,  y(t),  so  that 

Lt[x(t)|  =  y(t)  .  (2.1.1) 

ITic  operator  Ltj  ]  has  the  property  of  linearity;  that  is,  if 

Ljx,(t)]  =yi(t)  (2.1.2) 

and 

Ei|x2(t))  =  y2(t)  ,  (2.1.3) 

then 

Lilaxi(t)-t-Px2(t)l  =  ayi(t)  +  Py2(i)  (2.1.4) 

for  any  scalar  a  and  [).  However,  since  the  operator  Ld  |  evolves  over  time. 


Lt|x(t-t-X)|  ;<iy(t-t-x) 


(2.1.5) 


One  suitable  form  for  the  operator  Lt[  ]  is  the  discrete- time  superposition  sum 


y(t)  = 


oo 

X  k(t,x)  x(x)  , 


t  =  -  oo 


(2.1.6) 


where  the  function  k(t,x)  is  the  system  response  function,  or  Green’s  function  (see 
D'Angelo,  1970,  pp.  63-65).  In  one  sense,  the  system  response  function  characterizes  the 
response  of  the  system  at  some  time,  t,  to  an  input  applied  at  time  x.  The  system  can  also 
be  described  in  terms  of  a  time-varying  impulse  response,  h(t,A.),  as 


y(t)  = 


oo 

X  h(t,X)  x(t-X) 

X  ~  ~oo 


(2.1.7) 


The  time-varying  impulse  response  can  be  thought  of  as  the  response  of  the  system  at  time  t 
due  to  an  input  signal  applied  X  samples  earlier.  The  system  response  and  the  time-varying 
impulse  response  are  related  by 

k(t,x)  =  h(t,t-x)  .  (2.1.8) 

In  the  discussions  which  follow,  we  will  concern  ourselves  primarily  with  the 
system  response  function,  k(t,x).  This  choice  will  facilitate  our  presentation  of 
characterization  techniques  for  multirate  discrete-time  systems  (e.g.,  decimators  and 
interpolators).  In  these  systems,  the  input  and  output  sampling  rates  are  not  the  same.  Use 
of  the  time-varying  impulse  response  function.  h(t,X),  could  lead  to  ambiguities  in  the 
characterization  of  multirate  systems  since  it  is  not  clear  what  is  meant  by  "the  response  of 
the  system  at  time  t  due  to  an  input  signal  applied  X  samples  earlier." 

2.2  STABILITY  CONSIDERATIONS  FOR  LINEAR  TIME-VARYING 
SYSTEMS 

It  IS  relatively  straightforward  to  show  that  the  system  is  Bounded  Input  -  Bounded 
Output  (BIBO)  stable  for  each  t  if 


(2.2.1) 


4 


for  each  t.  First,  assume  that  k(t,t)  is  absolutely  summable  over  t  for  each  t  and  that  we 
apply  some  bounded  input,  x(t),  such  that 

lx(x)l<X<«>  (2.2.2) 

for  all  T.  It  then  follows  that 


ly(t,)l  = 


k(t,T)  X{T)I  < 


I 


Ik(t.x)  x(x)l 


(2.2.3) 


X  = 


I  = 


or 


ly(l)l  <  X 


^  Ikd.x 


)1< 


(2.2.4) 


X  =  -  OO 


Now  assume  k(t,x)  is  not  absolutely  .summable  over  all  x  for  some  i,  and  define 


x(x)  =  sgn  (k(t,x)) 


for  that  particular  t.  We  have  introduced  the  function 


(2.2.5) 


sgn(Y) 


+1  ;  Y  >  0 

- 1  ;  Y  <  0  ■ 


(2.2.6) 


Notice  that  x(x)  is  a  bounded  function  for  all  x.  The  system  output  corresponding  to  this 
input  is  given  by 


y(t)  = 


I 


k(t,x)  sgn(k(t.x))  = 


I 


lk(t,x)l 


X  =  -  -XJ 


X  = 


(2.2.7) 


which,  by  assumption,  is  unbounded. 

Hereafter,  we  will  assume  that  the  systems  of  interest  are  BIBO  stable  for  all  t;  in 
other  words 


s 


(2.2.8) 


holds  for  all  t  We  make  this  assumption  for  two  reasons.  First,  a  detailed  study  of  criteria 
for  the  stability  of  LTV  systems  is  beyond  the  scope  of  the  present  effort.  Second,  by 
making  this  reasonable  assumption  on  the  systems  of  interest,  we  obtain  interesting  and 
useful  results  when  we  consider  frequency-domain  techniques  for  characterizing  LTV 
systems.  For  a  more  detailed  discussion  of  stability  issues  for  LTV  systems,  see 
D’Angelo,  1970,  pp.  221-272. 

Notice  al.so  that  equation  (2.2.8)  implies 

ll<(t,T)l  =  0  (2.2.9) 

so  that  the  contribution  to  the  system  output  at  sample  time  t  due  to  an  input  at  some  sample 
time  t  in  the  remote  past  is  negligible.  In  a  sense,  then,  a  BIBO  stable  system  can  be  said 
to  possess  limited  memory  (see  Zadeh  and  Desoer,  1963,  pp.  45-46).  This  will  prove 
useful  to  us  when  we  consider  the  estimation  of  statistical  parameters  for  random  processes 
produced  by  stable  LTV  systems. 

2.3  FREQUENCY  DOMAIN  CHARACTERIZATIONS  OF  LINEAR  TIME- 
VARYING  SYSTEMS 

In  order  to  introduce  the  concept  of  a  frequency-domain  characterization  for  LTV 
systems,  let  us  first  consider  the  characterization  of  a  LTl  system  in  the  frequency  domain. 
For  the  LTI  system,  the  impulse  response  function  does  not  vary  with  time,  t,  but  is  a 
function  only  of  the  time  difference,  t-x  (see  Zadeh  and  Desoer,  1963,  p.  154).  The 
superposition  sum  for  the  LTI  system  is  therefore  given  by 


oo 

y(t)  =  X  h(t-T)  .\(T)  , 

t  -  -  oo 


(2.3.1) 


where  h(?i)  is  the  familiar  LTI  impul.se  response  function  (.see  Zadeh  and  Desoer,  1963, 
p.  156).  Assume  that  the  input  to  the  LTI  system  is  a  complex  sinusoid  of  frequency  a;  in 
other  words 


6 


x(x)  =  e 


lat 


(2.3.2) 


for  -71  <  a  <  7t.  Tlie  output  of  the  LTI  system  due  to  this  input  signaJ  is  given  by 


ya(0=  Zu  h(t-x)e 

T  =  -  oo 


lat 


(2.3.3) 


or  by  a  change  of  variables 


y^ct)  = 


I 

\  =  -< 


h(A.)  e 


-iaA. 


lai 


(2.3.4) 


Notice  that 

\hiX)  e'‘“^l  =  lh(>i)l  (2.3.5) 

for  all  values  of  a,  so  for  a  BIBO  stable  LTI  system  the  summation  enclosed  in  brackets  in 
equation  (2.3.5)  above  converges  uniformly  to  a  function  which  is  continuous  in  a  by  the 
Weierstrass  M-test  (see  Kaplan,  1962,  p.  169).  Let 


H(a)  = 


oo 


V  w-  X 

Zw  h(/,)  e 

X,  =  —  oo 


(2.3.6) 


so  that  the  response  of  the  LTI  system  due  to  a  complex  sinusoid  of  frequency  a  is  given 
by 


y^^(t)  =  H(a;  e'“‘  .  (2.3.7) 

rhc  tunction  H(a)  is  referred  to  as  the  frcqueney  response  function  of  the  LTI  system  (see 
Kaplan.  1962,  pp.  95-96).  Notice  that  in  .some  .sen.sc  the  frequency  response  function 
indicates  how  frequency  components  of  the  input  signal  are  mapped  into  frequency 
components  of  the  output  signal. 
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Now  consider  a  deterministic  input  signal,  x(t),  whose  Fourier  transform  exists 
and  assume  that  the  Fourier  transform  of  the  output  signal,  defined  by 


Y(co)  = 


Zu  y(0  e 


exists.  Substituting  for  y(t)  in  (2.3.8)  above  yields 


(2.3.8) 


Y(o))  = 


I  I 


h(t-t)  x(x)  e 


l  =  —  OOX  =  —  OO 


(2.3.9) 


or  by  a  change  of  variables,  and  assuming  the  system  is  BIBO  stable 


( 


Y(co)  = 


OK,* 

I 


x(x)  e 


t  =  -  oo 


Y  oo 


h(X)  e 


-icoX 


(2.3.10) 


so  that,  since  the  Fourier  transform  of  the  input,  x(x),  exists 


Y(o))  =  X(o))  H(a)) 


(2.3.11) 


(.see  Priestly,  1981,  p.  210). 

Equation  (2.3.11)  presents  one  of  the  fundamental  concepts  of  the  theory  of 
frequency-domain  characterizations  for  LTI  systems;  a  Fourier  transform  can  represent  a 
linear  convolution  of  series  in  the  time  domain  by  a  simple  product  of  functions  in  the 
frequency  domain.  For  this  rea.son,  a  Fourier  transform  is  said  to  be  a  compatible 
transform  for  the  class  of  LTI  systems  (see  D'Angelo,  1970,  pp.  275-276).  One  early 
approach  to  the  frequency-domain  characterization  of  LTV  systems  was  based  on  the  idea 
of  defining  new  transforms  which  would  be  compatible  transforms  over  some  restricted 
class  of  LTV  systems  (see  Aseltine,  1954,  and  Johnson  and  Kilmer,  1962).  For  example, 
the  Mellin  transfomi  defined  by 
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oo 


(2.3.12) 


is  a  compatible  transform  for  continuous  LTV  systems  governed  by  the  Euler-Cauchy 
differential  equation  (see  Aseltine,  1954).  One  obvious  shortcoming  to  this  approach  is 
that  specific  compatible  transforms  must  be  developed  for  each  specific  class  of  LTV 
system,  depending  on  the  time-varying  differential  equation  governing  that  class  of 
systems.  Also,  the  familiar  notion  of  "frequency"  can  be  lost  in  the  new  compatible 
transform  domain.  For  example,  it  is  not  obvious  how  the  parameter  |i  in  the  Mellin 
transform  of  equation  (2.3.12)  above  can  be  related  to  some  familiar  notion  of  "frequency" 
for  the  signal  y(t). 

An  alternative  approach,  introduced  in  Zadeh,  1950,  and  presented  for  continuous¬ 
time  LTV  systems  in  Zadeh,  1961,  Gersho,  1963,  and  Brikker,  1966,  and  presented  for 
discrete-time  LTV  systems  in  Claasen  and  Mecklenbrauker,  1982,  and  Crochiere  and 
Rabiner,  1983,  pp.  100-126,  is  based  on  the  idea  of  using  an  incompatible  transform  to 
obtain  a  characterization  of  the  LTV  system.  This  incompatible  transform  retains  the 
concept  of  "frequency"  of  the  input  and  output  signals.  In  order  to  introduce  this  concept, 
we  first  consider  applying  a  complex  sinusoidal  input  of  frequency  a  to  a  LTV  system.  Let 
the  response  of  the  LTV  system  due  to  this  sinusoidal  input  be  denoted  as  k(t,a).  By 
equation  (2.1.6),  k(t,a)  is  given  by 


K(t,a) 


k(t,x)  e 


lOT 


t  =  -  oo 


If  the  LTV  system  is  BIBO  stable  for  all  t,  then 


(2.3.13) 


<  oo 


T  =  -  oo 


for  all  t,  and  since 

lk(t,t)  e'”^l  =  lk(t,T)l 


(2.3.14) 


(2.3.15) 
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for  all  t,  by  the  Weierstrass  M-test  K(t,a)  exists  and  is  continuous  in  a  for  all  t  (see 
Kaplan,  1962,  p.  169).  Note  that  K(t,a)  is  simply  the  response  of  the  LTV  system  to  an 
applied  complex  sinusoid  of  frequency  a.  We  could  therefore  refer  to  K(t,a)  as  the  time- 
varying  frequency  response  function  for  the  LTV  system.  If  the  system  were  time- 
invariant,  we  know  from  equation  (2.3.7)  that  this  response  function  would  be  a  scaled 
version  of  the  original  complex  sinusoid.  In  Zadeh,  1950,  the  time- varying  system 
function,  H(t,a),  for  a  LTV  system  is  defined  to  be 

-ioti 

H(t,a)  =  K(t,a)  e  (2.3.16) 

In  a  sense,  the  time-varying  system  function  is  the  complex  envelope  of  the  response  of  the 
LTV  system  to  an  applied  complex  sinusoid  (see  Zadeh,  1950). 

Now  consider  a  signal,  x(x),  whose  Fourier  transfomi  exists;  in  other  words, 

oo 

X(a)  =  ^  x(T)  e (2.3.17) 

X  =  -  oo 

exists.  We  could  therefore  describe  the  signal  x(t)  as 

jt 

x(t)-^  I  X(a)  e'“^  da  .  (2.3.18) 

-n 

Applying  this  signal  to  the  LTV  system  of  equation  (2.1.6)  yields  an  output  signal,  y(t), 
given  by 


yw  =  2K  ^ 

X  =  -  oo  -71 


lat 

e  da  . 


Interchanging  the  order  of  summation  and  integration  yields 


■n 


)  e 


lai 


X  =  -  oo 


X(a)  da 


(2.3.19) 


(2.3.20) 
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Notice  that  the  term  in  brackets  in  equation  (2.3.20)  above  is  the  time-varying  frequency 
response  function,  K(t,a).  Equation  (2.3.20)  may  therefore  be  rewritten  as 

n 

K(t,a)  X(a)  da  .  (2.3.21) 

Now  assume  that  the  Fourier  transform  of  the  output  signal,  y(t),  exists;  in  other  words 


=  I 


-iwt 


Y(oj)  =  ^  y(t)  e 

t  =  —  oo 


exists.  Substituting  equation  (2.3.21)  into  (2.3.22)  above  yields 


(2.3.22) 


Y({n)  = 


2k 


X  J  K(t, 


a)  X(a)  da  e 


10)1 


(2.3.23) 


I  s:  —  oo  ~1t 


Interchanging  the  order  of  summation  and  integration,  we  obtain 


-K 


K(t,a)  e 


-10)1 


X(a)  da 


(2.3.24) 


For  the  moment,  assume  that  the  summation  in  brackets  in  equation  (2.3.24)  above  exists, 
and  define 


K(to.a)=  x'(t,a)  e'‘“‘  (2.3.25) 

i  ~  oo 

Conditions  for  the  existence  of  K(ca,a)  and  the  implications  of  these  conditions  will  be 
examined  in  Section  2.5.  Substituting  equation  (2.3.25)  into  equation  (2.3.24)  we  obtain 


K 


K(a),a)  X(a)  da 


(2.3.26) 


II 


Notice  that  equation  (2.3.26)  in  some  sense  gives  a  description  of  how  input  frequency 
components  contribute  to  output  frequency  components.  It  is  not  the  simple  one-to-one 
relationship  one  finds  for  LTI  systems,  and  it  is  for  this  reason  that  this  characterization  of 
LTV  systems  is  considered  to  be  an  incompatible  transform  technique  (see  D'Angelo, 
1970,  p.  294).  The  function  K((o,a)  is  referred  to  as  the  time-varying  system 
transmission  function  (see  Claasen  and  Mecklenbrauker,  1982,  and  Crochiere  and  Rabiner, 
1983,  p.  103). 

2.4  ANALYSIS  OF  NETWORK  STRUCTURES  FOR  DISCRETE-TIME 
LINEAR  TIME-VARYING  SYSTEMS 


We  have  introduced  time  and  frequency-domain  representations  for  discrete-time 
LTV  systems  in  the  previous  sections.  In  this  section,  we  consider  the  analysis  of  network 
structures  obtained  by  connecting  LTV  .subsystems.  Two  fundamental  structures  will  be 
considered:  the  parallel  connection  of  two  LTV  subsystems,  and  the  series  connection  of 
two  LTV  subsystems.  More  complicated  network  structures  may  be  analyzed  by  extension 
of  the  methods  presented  below. 


Consider  two  LTV  systems  with  response  functions  ki(t,x)  and  k2(t,x).  Assume 
that  the  corresponding  transmission  functions,  Ki((0,a)  and  K2(ti),a),  exist  for  these  two 
systems.  A  parallel  connection  of  the  systems  is  indicated  schematically  in  Figure  2.4.1. 
The  response  of  the  overall  .system  to  some  input  signal,  x(t),  is 


y(t)  = 

X  ki(t,T)x(T) 

+ 

X  k2(t,T)x(T) 

X  -  -  OQ 

K  ) 

(2.4.1) 


If  we  were  to  apply  a  unit  impulse  function  at  some  time  T  =  Tq  to  this  composite  system, 
the  observed  output  would  be  the  response  function  for  the  overall  LTV  system  for  the 
fixed  input  time  i  =  To,  k(t,To).  l.ot  the  unit  impul.se  function  be  denoted  by  ^(t);  in  other 
words  define 


C(T)  = 


t  =  0 
T  0 


(2.4.2) 


The  response  function  for  the  overall  system  is  given  by 
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k(t,to)  = 


f  <30  \ 

s 


ki(t,t)^(T-to) 


X  =  -oo 


k2(t,X)C('C-to) 


t  =  —  oo 


Simplifying  this  expression  we  obtain 


k(t,To)  =  kutJo)  +  k2(t,To)  ; 


(2.4.3) 


(2.4.4) 


in  other  words,  the  response  function  of  a  parallel  connection  of  LTV  systems  is  simply  the 
sum  of  the  individual  response  functions.  Also  note  that  if  ki(t,x)  and  k2(t,t)  represent 
BIBO  stable  LTV  systems,  then  the  overall  LTV  system  is  also  BIBO  stable. 


Next  consider  applying  a  complex  sinusoid  ot  frequency  a  to  the  composite  system 
of  Figure  2.4.1.  We  know  that  the  output  of  the  composite  system  for  this  input  will  be 
the  frequency  response  function  for  the  overall  system  at  frequency  a,  K(t,a).  From 
equation  (2.4. 1 )  we  may  express  K‘(t,a)  as 


K(t,a) 


(t,t)  e 


KIT 


k2(t,x)  e 


lOT 


t  =  -  oo 


(2.4.5) 


Simplifying  this  expression,  we  obtain 


hJtt.a)  =  Ki(t,a)  +  K2(t,a)  , 


(2.4.6) 


where  x'ltt.a)  and  K2(t,(x)  arc  the  frequency  response  functions  corresponding  to  ki(t,x) 
and  k2(t,X).  From  equation  (2.3.25)  we  can  express  the  transmission  function  of  the 
overall  system,  K((o,a),  as 


K( 


I 


(0,u)  =  (Ki(t,a)  +  x'2(t,a))  e 


10)1 


(2.4.7) 


Simplifying  this  expression  we  obtain 


Kt(o.«)  =  K|((o.(t)  +  K2(to.(x)  , 


(2.4.8) 
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where  Ki(ca,a)  and  K2(o),a)  are  the  transmission  functions  corresponding  to  ki(t,T)  and 
k2(t,t).  From  equation  (2  4.8)  we  see  that  the  transmission  function  for  a  parallel 
connection  of  LTV  systems  is  simply  the  sum  of  the  individual  transmission  functions. 


Finally,  consider  the  series  connection  of  LTV  systems  indicated  schematically  in 
Figure  2.4.2.  The  response  of  the  overall  system  to  .some  input  signal,  x(t),  is  given  by 


=  I 


y(t)=  ^  k2(t,Y)yi(Y) 

Y  =  -  « 


(2.4.9) 


where 


yi(Y) 


=  X  ki(Y.t) 


(2.4.10) 


Combining  equations  (2.4.9)  and  (2.4.10),  we  obtain 


y(t)=  k2(t,Y)  ki(Y,x)x(x)  , 

Y=-oo  T  =  -'*> 


(2.4.11) 


or  simply 


y(t)=  Z  Z 

t  =  -  oo  I  Y  =  -  t 


k2(t.Y)  l<i(TX)  x(x) 


(2.4.12) 


If  we  apply  a  unit  impulse  at  time  x  =  Xo  to  this  composite  system,  we  obtain  the  overall 
system  response  function  for  input  time  x  =  Xq  .  k(t,Xo).  From  equation  (2.4.12),  we 
may  express  k(t,Xo)  as 


Z  Z 


k2(t.Y)  ki(TX)  k(^-'to)  . 


(2.4.13) 
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or  simply 


t.to)  =  X 


k2(t.Y)  ki(Y,To) 


(2.4.14) 


Equation  (2.4.14)  indicates  that  the  overall  system  response  function  at  input  time  t  =  Xq 
for  a  series  connection  of  two  LTV  systems  is  the  response  of  the  second  LTV  system 
when  excited  by  the  first  LTV  system's  response  function,  Ici(y,Xo)-  Notice  that,  in 
general,  a  different  response  function  (and  therefore  a  different  LTV  system)  is  obtained 
when  we  commute  the  systems  kt(t,x)  and  k2(t,T)  (see  Crochiere  and  Rabiner,  1983, 
pp.  108-1 12).  Also  note  that  if  both  ki(t,t)  and  k2(t,x)  are  BIBO  stable  LTV  systems, 
then  the  overall  system  is  also  BIBO  stable. 

Now  consider  applying  a  complex  exponential  of  frequency  a  to  the  composite 
LTV  system  of  Figure  2.4.2.  The  output  signal  obtained  in  this  case  will  be  the  frequency 
response  function  for  the  overall  system  at  frequency  a,  K(t,a).  From  equation  (2.4.12) 
we  may  express  K(t,a)  as 


.,a)=  I  I 


T  =  -  oo  Y  =  - 


k2(t,Y)  ki(Y.t)  e 


(2.4.15) 


On  simplifying  equation  (2.4.15),  we  obtain 


t,a)  =  ^  k2{t,Y)  Ki(Y,a)  . 


(2.4.16) 


Consider  that  we  may  express  k2(t,Y)  as 


k2(t,Y)  =  2^  1  K'2(t.(J)  e  dp  , 


(2.4.17) 


so  that  equation  (2.4.16)  can  be  rewritten  as 
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(2.4.18) 


n 

K(t,a)  =  ^  J  K2(t,p)  X  Ki(Y,a)  e 

"K  'y  s:  —  oo 


dp 


Notice  that  the  expression  in  brackets  in  equation  (2.4.18)  above  is  the  transmission 
function  for  the  first  LTV  system,  Ki((0,a).  Our  expression  for  K(t,a)  therefore  becomes 


K:(t,a)  =  ^  j  K:2(t,p)  Ki(p,a)  dp  . 


(2.4.19) 


Using  equation  (2.3.25)  to  obtain  the  transmission  function  for  the  overall  LTV  system, 
K(co,a),  we  obtain 


K((o,a)  = 


t  =  -  oo 


^  f  ^ 

X  ^  j  »C2(t,p)  Ki(P,a)  dP 

V 


-jt 


-10)1 


(2.4.20) 


or  simply 

K 

K(co,a)  =  ^  j  K2(co,P)  Ki(P,a)  dp  .  (2.4.21) 

-71 

Recall  that,  in  some  sense,  Ki(p,a)  indicates  the  distribution  of  frequency  components  in 
the  output  of  the  first  LTV  system  due  to  a  complex  exponential  input  signal  of  frequency 
a.  Comparing  equation  (2.4.21)  and  (2.3.26),  we  see  that  the  overall  system's 
transmission  function,  K(ti),a).  is  simply  the  response  of  the  second  LTV  system  due  to  a 
signal  with  frequency  content  Ki(P,a).  Again,  note  that  in  general  commuting  the  two 
LTV  systems  yields  an  overall  LTV  system  with  a  different  transmission  function  (see 
Crochiere  and  Rabiner,  1983,  p.  111). 
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FIGURE  2.4.1 

PARALLEL  CONNECTION  OF  LTV  SYSTEMS 


AS-9 1-364 


AS-91-365 


FIGURE  2.4.2 

SERIES  CONNECTION  OF  LTV  SYSTEMS 
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2.5  CLASSES  OF  DISCRETE-TIME  LINEAR  TIME-VARYING  SYSTEMS 


The  discussions  of  frequency-domain  characterizations  of  LTV  systems  using  the 
system  transmission  function,  K(co,a),  presented  thus  far  have  been  made  by  simply 
assuming  that  K((o,a)  exists.  We  now  turn  our  attention  to  the  conditions  under  which 
K(a),a)  exists  and  the  implications  of  these  conditions  on  the  properties  of  the  LTV 
system.  Recall  that  the  system  transmission  function  was  defined  to  be 


-icot 

K(a},a)  2^  K(t,a)  e  (2.5.1) 

t  =  -  oo 

Clearly,  the  existence  of  K(OJ,a)  intimately  depends  on  the  properties  of  the  frequency 
response  function,  K(t,a).  We  will  consider  two  classes  of  LTV  system:  the  class  of  LTV 
systems  whose  frequency  response  function  is  a  periodic  function  of  t  for  all  a,  and  the 
class  of  LTV  systems  whose  frequency  response  function  is  an  almost  periodic  function  of 
t  for  all  a. 

First  assume  that  the  frequency  response  function  is  a  periodic  function  of  output 
time  t  for  all  a;  in  other  words,  assume 

K(t,a)  =  K(t+T(a),a)  (2.5.2) 

for  some  T(a)  for  all  a.  Notice  that  equation  (2.5.2)  assumes  that  the  frequency  response 
function  may  be  periodic  with  a  distinct  period,  T(a),  for  each  a.  Since  K(t,a)  is  a 
periodic  function  in  t,  it  has  a  Fourier  series  representation  given  by 

oo 

Zi27:nt/r((x) 

an(a)  e  (2.5.3) 


w  here 


an(a) 


'1  KZ) 

I 


1=0 


i2nniyT(n) 

K'd.fX)  e 


(2.5.4) 
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(see  Oppenhiem  and  Schafer,  1975,  pp.  88-89).  Substituting  equation  (2.5.3)  into 
(2.5.1),  we  obtain 


K(co,a)  = 


S  I 


i27tnt/T(a>  -icot 

an(a)  e  e 


=  —  CX5  —  OO 


(2.5.5) 


or  simply 


K(a) 


V  ,  ,  V  -i(a)-2jtn/T(a))l 

,a)  =  ^  an(a)  ^  e 


t  =  -  oo 


(2.5.6) 


The  temi  in  brackets  in  equation  (2.5.6)  is  given  by 


Ic 


t  =  -  oo 


-i(a)-27ui/^r(a)H 


=  6(to-2jtn/T(tt)) 


(2.5.7) 


where  0(7)  is  the  Dirac  delta  function  (see  Papoulis,  1962,  pp.  42-46).  Our  expression 
for  Ktto.cx)  is  therefore  given  by 


K(a),a) 


I 


an((x)  5(a)-27tnAr(a))  . 


n  =  -  00 


(2.5.8) 


From  equation  (2.5.8)  above  w'c  see  that  a  Ll'V  system  whose  frequency  response 
function  is  a  periodic  function  of  time  t  will  possess  a  system  transmission  function  which 
consists  of  niixlulated  impulse  sheets  in  the  (w.a)  bifrequency  plane. 

Finally,  assume  that  the  frequency  re.spon.se  function  is  an  almost  periodic  function 
of  output  time  t  for  all  tt;  in  other  words,  for  any  e  >  0,  there  is  some  T(a)  such  that 

Ixit.tx)  -  k(  t-t-T((ii.(x)l  <  r  (2.5.9) 
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for  all  a  (see  Wiener,  1933,  p.  185).  Since  K(t,a)  is  an  almost  periodic  function  for  all  a, 
we  know  ihat  there  is  an  increasing  sequence  of  frequencies, 
{.  .  .,  a)-i(a),  coo(a),  toi(a), .  .  such  that 


oo 

K(t,a)  =  51  an(a) 
n=  — oo 


(2.5.10) 


where 


an(a)  = 


lim  j_ 
to~^°°  2to 


K:(t,a) 


-ia)_(a)i 

e 


(2.5.11) 


(see  Wiener,  1933,  p.  186  and  Priestly,  1981,  p.  199).  Substituting  equation  (2.5.10) 
into  (2.5.1),  we  obtain 


K((o,a) 


Z  Z  a„(a>  e  •“ 


(2.5.12) 


or  simply 


K(o),a)  =  2^  an(a)  6(a)-c0n(a))  .  (2.5.13) 

n  =  — oo 

Again,  we  see  that  LTV  systems  which  possess  almost  periodic  frequency  response 
functions  exhibit  transmission  functions  consisting  of  modulated  impulse  functions  in  the 
(a),a)  bifrequency  plane. 

In  Claa.sen  and  Mecklenbrauker,  1982,  the  class  of  LTV  systems  with  almost 
periodic  frequency  response  functions  is  referred  to  as  the  class  of  "stationary"  LTV 
systems.  The  authors  then  point  out  that,  whereas  series  and  parallel  connections  of  LTV 
systems  with  periodic  frequency  response  functions  may  not  yield  overall  LTV  systems 
with  periodic  frequency  response  functions,  series  and  parallel  connections  of  LTV 
systems  with  almost  periodic  frequency  response  functions  always  yield  overall  LTV 
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systems  with  almost  periodic  frequency  response  functions.  In  the  discussions  which 
follow,  we  will  refer  to  the  class  of  LTV  systems  with  almost  periodic  frequency  response 
functions  as  simply  the  class  of  linear  periodically  time  varying  (LPTV)  systems. 

2.6  IMPLICATIONS  OF  CAUSALITY  ON  THE  TRANSMISSION 
FUNCTION  OF  LINEAR  PERIODICALLY  TIME-VARYING 
SYSTEMS 


Recall  from  equation  (2.5.13)  that  the  system  transmission  function  for  a  LPTV 
system  is  given  by 


K((o,a) 


I 


an(a)  5((o-con(a)) 


n  =  - 


(2.6.1) 


The  system  transmission  function  is  therefore  non-zero  only  along  the  set  of  curves 
O)=C0n(a)  in  the  (a),a)  bifrequency  plane.  In  Claasen  and  Mecklenbrauker,  1982,  it  is 
shown  that  for  causal  continuous-time  LPTV  systems,  the  system  transmission  function  is 
non-zero  only  along  a  set  of  lines  in  the  bifrequency  plane.  It  is  then  asserted  that  causal 
discrete-time  LPTV  systems  will  similarly  possess  O)n(a)  which  are  linear  in  a.  In  this 
section,  we  will  show  that  causal  LPTV  systems  whose  input  and  output  sampling  rates  are 
rationally  related  do  indeed  have  the  property  that  the  function  a)n(a)  is  linear  in  a. 

First,  consider  a  one-sided  discrete-time  signal,  s(y);  in  other  words,  a  sequence 
for  y  =  ...  ,  -1,  0,  1,  ...  such  that 

s(y)  =  0,  Y<0  .  (2.6.2) 

Assume  that  the  Fourier  transfomi  of  s(y)  exists  and  is  given  by 


oo 


Y  =  0 


(2.6.3) 


Now  consider  the  real  and  imaginary  components  of  S(Q),  given  by 
S(Q)  =  Sr(Q)  +  iSi(i2)  . 


(2.6.4) 
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For  a  one-sided  signal,  the  real  and  imaginary  components  of  S(Q)  are  related  through  a 
discrete  Hilbert  transform  (see  Oppenheim  and  Schafer,  1975,  pp.  339-345).  Specifically, 
the  real  and  imaginary  components  of  S(Q)  are  related  by 

n 

Si(Q)  =  ^(p.v.)  I  SR(0)cot^^jd0  (2.6.5) 

and 

rt 

SR(a)  =  -  ^(p.v.)  j  Si(0)cot|^^jd0  +  s(O)  .  (2.6.6) 

In  equations  (2.6.5)  and  (2.6.6),  (p.v.)  indicates  the  Cauchy  principle  value  of  the 
integral.  In  one  sense,  the  Hilbert  transform  is  a  circular  convolution  of  two  functions  in 
Q,  and  so  is  itself  a  linear  function  in  Q  (see  Oppenheim  and  Schafer,  1975,  p.  345). 

We  must  now  define  precisely  what  is  meant  by  a  causal  LTV  system.  First, 
assume  that  the  input  and  output  sampling  times,  T  and  t,  can  be  related  to  some  underlying 
continuous  time  parameter,  r.  Specifically,  we  will  assume  that  output  samples  are  taken  at 
unit  intervals  in  /,  while  input  samples  are  obtained  at  intervals  spaced  p  apart.  This  is 
indicated  schematically  in  Figure  2.6.1.  By  coupling  input  and  output  sampling  times  to 
this  underlying  continuous  parameter,  we  can  make  some  sensible  accounting  of  which 
output  samples  occur  after  a  given  input  sample.  A  LTV  system  with  system  response 
function  k(t,t)  is  defined  to  be  causal  if 

k(t,x)=0,  t<[pT)  .  (2.6.7) 

The  notation  (x)  for  some  real-valued  x  indicates  the  smallest  integer  which  is  greater  than 
or  equal  to  x.  Notice  that  equation  (2.6.8)  indicates  that  ob.served  outputs  at  time  t  depend 
only  on  inputs  applied  at  some  earlier  absolute  time.  Systems  for  which  p  =  1  will  be 
called  constant-rate  systems,  since  the  input  and  output  sampling  rates  are  the  same. 
Systems  for  which  p  <  1  or  p  >  1  will  be  called  rate-reduction  and  rate-increase  systems, 
respectively. 
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Assume  that  k(t,T)  is  the  system  response  function  for  a  causal,  stable,  constant- 
rate  LPTV  system.  From  equations  (2.3.13)  and  (2.3.25),  we  may  express  the  system 
transmission  function  as 


oo  oo 


).«)=  X  k(t. 


,  lat  -Kot 
X)  e  e 


1=  —  OOX  =  — oo 


Now  consider  that 


,  V  w  ,  i(a)-a)T  -icot 

K(to,to-a)  =  Za  k(t,t)  e  e 


l  =  — oo  X  = 


(2.6.8) 


(2.6.9) 


oo  /  oo 


K(o),(o-a) 


=  I  z 


X  =  — oolts;— oo 


k(.,T)  e-'“' 


(2.6.10) 


Substituting  Y=  t-t,  we  obtain 


K(co,a)-a)  =  ^  ^  k(y+t,t 

X”— oo  y=— oo 


^  -icoy  -lax 
)  e  e 


(2.6.11) 


Define  the  sequence  s(Y,t)  to  be 


s(Y,T)  =  k(Y+t,t)  , 


(2.6.12) 


so  that  equation  (2.6.1 1)  becomes 


K(a),a)-a)  =  ^  ^ . 

t  =  -  oo  Y  =  —  ' 


-iwy  -lat 
s(Y,t)  e  e 


(2.6.13) 


23 


The  summation  in  brackets  in  equation  (2.6.13)  above  is  the  Fourier  transform  of  s(y,x) 
with  respect  to  y.  However,  from  equation  (2.6.10)  we  know  that  s(y,T)  is  a  one-sided 
sequence;  in  other  words, 

s(y,t)  =  0  ,  y  <  0  .  (2.6.14) 

If  we  define 


S(co,x) 


Y=  -oo 


=  SR(a),x)  +  iSi(a),x)  , 


(2.6.15) 


then  we  know  that  the  real  and  imaginary  components  of  S((j0,x)  are  related  by 


Si(to,x)  =  ~  (p.v.)  J  Sr(0,x)  cot 


^0)-© 


-t: 


d© 


and 


SR(to.x)  =  - 


^  (p.v.)  J  Si(©,x)  cot 


d©  +  k(x,x)  . 


(2.6.16) 


(2.6.17) 


Substituting  equation  (2.6.15)  into  equation  (2.6.13)  and  simplifying,  we  obtain 


K(ct),a)-a)  =  2^  Sr(o),x)  cos(ax)  -i-  Si((o,x)  sin(ax) 

X  =  -  oo 


■  I 


X  =  -  oo 


Si(co,x)  cos(ax)  -  Sr((o,x)  sin(ax) 


(2.6.18) 


From  equation  (2.6.18),  we  may  write  the  real  and  imaginary  components  of  the 
transmission  function  as 


KR(o),aj-a)  =  SrCo),!)  cos(ax)  +  Si(a),x)  sin(ax) 

t  =  -  oo 


(2.6.19) 
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and 


=  I 


Ki(o),a)-a)  =  ^  Si(o),t)  cos(at)  -  SR(a),x)  sin(aT)  , 

X  =  -  oo 


(2.6.20) 


respectively.  Now  consider  that 


^(p.v.)  f  KR(0,0-a)cot[^^jd0  = 

^  ^(p  v.)  J  SR(0,t)  cot  ^^ld0cos(aT) 

X  =  —  oa  -K  \  J 

+  S  ^  (p.v.)  /  Si(0,t)  cot  ® 


X  =  -  oo 


-7t 


V  ^  y 


d0  sin(ax) 


(2.6.21) 


or,  by  equations  (2.6.16)  and  (2.6.17), 


It 

^(p.v.)  j  KR(0,0-a)cotj^^jd0  = 


2^  Si(to,x)  cos(ax)  -  SR(a),x)  sin(ax)  +  Z,  k(x,x)  sin(ax)  .(2.6.22) 

x  =  — oo  X  =  — 


From  equations  (2.6.22)  and  (2.6.20),  we  obtain 


Ki(a),o)-a)  =  (p.v.)  J  KR(0,0-a)cotf^ld0-Ho(a),  (2.6.23) 

-It  \  ^  J 


where 


Ho(a)  =  2^  k(x,x)  sin(ax) 

X  =  -  OO 


(2.6.24) 


an  odd  function  of  a.  Proceeding  in  a  similar  fashion,  we  detemiine  that 
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(2.6.25) 


71 

KR(o),Q)-a)  =  -  ^  (p.v.)  /  Ki(e,0-a)  cot  +  He(a) , 


where 


He(a)  = 


X 


k(T,T)  cos(ax)  , 


X  =  -  oo 


(2.6.26) 


an  even  function  of  a.  From  equation  (2.6.8),  it  is  straightforward  to  show  that 

7t 

J  K(co,w-a)  dco  =  He(a)  -  iHo(a)  .  (2.6.27) 

-n 

It  should  also  be  noted  that  equations  (2.6.23)  through  (2.6.27)  hold  for  any  causal, 
constant-rate  LTV  system  for  which  K(a),a)  exists. 

Consider  equation  (2.6.1)  and  note  that  we  may  express  the  real  and  imaginary 
components  of  the  transmission  function  as 


KR(o),o)-a) 


I 


aRn((o-a)  5(<i)  -  (Dn(o)-a)) 


n=  -oo 


(2.6.28) 


and 


Ki((o,a)-a) 


X 


a[n((o-a)  5(to  -  tOn(o)-a))  , 


n=  -oo 


(2.6.29) 


where  aRn(oc)  and  ain(a)  are  the  real  and  imaginary  components  of  an(cx).  From 
equation  (2.6.23),  we  obtain 
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ain(o)-a)  d(o)  -  a)n(o>-a))  = 


“  n 

S  ^(Pv)l  aR„(e-a)  5(0-cOn(0-a))  cot  f^ld0-Ho(a).  (2.6.30) 

n  —  —  oo  V  / 


Notice  that  the  left-hand  side  of  equation  (2.6.30)  is  identically  zero  for  all  O)  ^  a)n(to-a), 
which  implies  that  the  summation  on  the  right-hand  side  of  equation  (2.6.30)  is  identically 
zero  for  all  to  ^  a)n(to-a).  This  will  occur  only  if  5(o)-C0n((o-a))  is  not  a  function  of  to, 
in  other  words  only  if  tOn(a)  is  a  linear  function  of  a.  We  can  therefore  express  a)n(a)  as 


con(a)  =  a  4  Won  . 


(2.6.31) 


where  Won  is  some  initial  offset  for  the  n-th  response  curve.  From  equations  (2.6.1), 
(2.6.27),  and  (2.6.31),  we  obtain  the  relation 


^  J  K(w,w-a)  dw  =  An  5(a-Won)  . 


(2.6.32) 


where 


An  -  J  3n(P)dP 


(2.6.33) 


so  that 


Hc(a)  =  2^  Aro  6(a-Won) 


(2.6.34) 


Ho(a)  = 


A|n  S(a-Won) 


(2.6.35) 
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where  Arh  and  Arh  are  the  real  and  imaginary  components  of  An.  To  summarize,  we  have 
seen  that  a  causal,  constant-rate  LTV  system  with  an  almost  periodic  frequency  response 
function  has  a  system  transmission  function  given  by 


K(a),a)  =  an(a)  5(o)  -  a  -  tOon)  •  (2.6.36) 

n  =  -oo 

From  equation  (2.6.36)  above  we  see  that  the  transmission  function  for  a  causal,  constant- 
rate  LPTV  system  is  non-zero  only  along  a  set  of  lines  with  unit  slope  in  the  (0),a)  plane. 
The  set  of  offset  frequencies,  (Oon.  in  one  sense  can  be  interpreted  as  the  response  of  the 
LPTV  system  due  to  a  constant  unit  input  sequence  and  is  therefore  determined  by  K(t,0). 

Now  consider  the  case  of  a  causal,  integer  rate-increase  LPTV  system,  in  other 
words  a  system  for  which  p  =  m,  an  integer.  From  equation  (2.6.8),  we  have  the 
relation 


oo  oo 

V,  X  ^  ^  /  X  i(m(i)-a)-c  -iwi 

K(a),ma)-a)  =  Zj  Zw  k(t,x)  e  e 

l=—  OOXS—  oo 


(2.6.37) 


or  by  a  change  of  variables 


K(a),m(o-a) 


I 


f 


I 


s(Y.'t) 


-itOY 

e 


-iUT 

e 


X  =  -  oo 


Y  =  -' 


J 


In  equation  (2.6.38),  we  have  defined  the  one-sided  sequence 


(2.6.38) 


s(Y,T)  =  k(Y+mx.t)  =  0  ,  Y  <  0  .  (2.6.39) 

The  term  in  brackets  in  equation  (2.6.38)  is  the  Fourier  transform  of  s(Y,x)  with  respect  to 
Y,  so  that  the  Hilbert  transfonn  relations  of  equations  (2.6.16)  and  (2.6.17)  still  apply. 
Proceeding  along  the  same  lines  presented  for  the  constant-rate  system,  we  obtain 
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7t 

Ki(to.mci>-a)  =  ~  (p.v.)  J  KR(0,ni0-a)  cot  *  Ho(a) ,  (2.6.40) 


where 


Ho(a)  =  Zj  k(mx,t)  sin(ax)  , 


(2.6.41) 


KR(0),ma)-a)  =  -  ^  (p.v.)  J  Ki(©,m©-a)  cot  ^^^)d©  +  He(a) ,  (2.6.42) 


where 


He(a)  =  ^  k(mx,x)  cos(ax) 


(2.6.43) 


Furthermore,  it  is  straightforward  to  show  that 


J  KCco.mco-a)  dco  =  Hc(a)  -  iHo(a) 


(2.6.44) 


It  is  appropriate  to  note  that  equations  (2.6.40)  through  (2.6.44)  apply  for  any  causal, 
integer-rate-increase  LTV  system  for  which  K(o),a)  exists. 

From  equation  (2.6,1),  we  may  write  the  real  and  imaginary  components  of 
K(o),m^o-^a)  as 


KR(to,mto-a) 


~  aRn 


(mo)-a)  5(co  -  ton(mto-a)) 


(2.6.45) 


Ki(to,mtD-a)  = 


(2.6.46) 


ain(ma>-a)  5(co  -  tOn(mo)-a))  , 


where  aRn(ct)  and  ain(a)  are  the  real  and  imaginary  components  of  an(a).  From 
equation  (2.6.40),  we  obtain 


ain(ma)-a)  5(a)  -  con(ma)-a))  = 


^  ^p  v  )  j  aRn(m0-a)5(0-O)n(m©-a))cot  p0  -  Ho(a) .  (2.6.47) 

n  =  -  oo  -n  V  / 


We  can  again  argue  that  the  only  way  that  both  sides  of  equation  (2.6.47)  above  can  be 
zero  for  o)  ^  O)n(ma)-a)  is  for  5(co  -  o)n(mo)-a))  to  be  independent  of  O).  This  is  the 
case  only  if  a)n(a)  is  a  lire  j  function  such  that 


Wn(a)  -  In  +  0)on  , 


(2.6.48) 


wher<"  0)on  is  some  initial  offset  for  the  n-th  response  curve.  From  equations  (2.6.1), 
(2.0.44),  and  (2.6.48)  we  obtain 


K  ^ 

~  J  K(a),ma)-a)  do)  =  2^  An  5(a-a)on)  , 


(2.6.49) 


where 


j  an((3)  dp 


(2.6.50) 


so  that 


Hc(a)  = 


Aro  5((x-a>^)n) 


(2.6.51) 


.^0 


and 


Ho(cx)  —  ^  Ajn  5(a-ci)on)  >  (2.6.52) 

n  =  -  oo 

where  Arh  and  Ain  are  the  real  and  imaginary  components  of  An-  To  summarize,  we  have 
seen  that  a  causal,  integer-rate-increase  LPTV  system  has  a  system  transmission  function 
given  by 


K(o),a) 


(2.6.53) 


From  equation  (2.6.53)  above  we  see  that  the  transmission  function  for  a  causal,  integer- 
rate-increase  LPTV  system  is  non-zero  only  along  a  set  of  lines  with  slopie  1/m  in  the  (a),a) 
plane.  Again,  the  set  of  offset  frequencies,  Won.  may  be  determined  from  K(t,0),  the 
response  of  the  system  to  a  constant  unit  input  sequence. 


Finally,  assume  that  k(t,T)  represents  a  causal,  stable  LPTV  system  whose  input 
and  output  sampling  rates  are  rationally  related,  in  other  words  p  =  m/n  for  some  integer  m 
and  n.  For  the  moment,  consider  that  we  could  rewrite  equation  (2.6.8)  as 

n- 1  °°  oo 

K(co,a)  =  ^  ^  k(t,ny-Hj)  e (2.6.54) 

j  =  0  y=-ool=:-oo 

Now,  for  j  =  0,  1,  •  ,  n- 1.  define 


tj  = 


njj 

n 


so  that,  for  P  =  t  -  tj  ,  equation  (2.6.54)  can  be  rewritten  as 


(2.6.55) 


V’  ^  ,  r,  '“(ny+j)  -iw(P+tj) 

ZZ  kiP-t-t,,  ny+i)  e  e 


n-l 

K(o).a)  =  I  Z^  ZLu  k(P-t-tj,  ny+j)  e 

j  =  0y  =  -oo(3  = 


(2.6.56) 


Rearranging  terms  in  equation  (2.6.56),  we  obtain 


n-1  «  S  ^  .o  V  ‘"“Y  ‘“j  -‘“‘j 

K(o),a)  =  £  S  S  Kj(P,Y)  e  e  e  e 
j  =  OlY=-“P  =  -“’ 


(2.6.57) 


where  we  have  defined  the  function 


lcj(p,Y)  =  k(b+tj,  nY+j)  . 


(2.6.58) 


A  great  deal  of  insight  into  the  behavior  of  the  system  represented  by  ^j(P,Y)  can  be 
obtained  by  a  careful  scrutiny  of  equation  (2.6.58).  First,  consider  that  for  each  j,  fcj(P,Y) 
is  the  system  response  function  for  a  stable  LPTV  system.  Stability  is  obtained  from  the 


fact  that 


(2.6.59) 


is  a  subseries  of 


(2.6.60) 


which  converges  for  each  p.  Similarly,  it  is  straightforward  to  show  that,  since 


I  =  -lo  !=:-<« 


(2.6.61) 


exists  for  some  set  of  tOn(a), 


P  =  -l«  Y  =  -  ' 


(2.6.62) 


also  exists  for  the  same  set  of  o>n(a).  Finally,  consider  that 
fcjlp.y)  =  0  for  all  P  <  my  ; 

in  other  words.  (sjiP.y)  is  a  causal,  integer-rate-increase  system. 


(2.6.63) 


The  double  summation  in  equation  (2.6.57)  above  is  recognized  as  the  system 
transmission  function  corresponding  to  (tj(P,Y).  denoted  Kj(w,na).  We  therefore  obtain 


(2.6.64) 


n-l 

Z>'  iaj  -icotj 

Kj(0),na)  e  e 

J  =  0 

From  the  discussion  of  in tege. -rate-increase  LPTV  systems  presented  above,  we  know  that 

A 

Kj(a),0)  is  non-zero  only  along  a  set  of  lines  defined  by 

(0  =  ®  -t-  Won  •  (2.6.65) 

From  equation  (2.6.62),  we  know  that  the  cOon  are  the  same  for  each  integer-rate-increase 
system  since  Won  =  Wn(fi)-  By  inspection  of  equation  (2.6.64),  we  may  therefore 
conclude  that  the  system  transmission  function  for  the  m/n  rational-rate  system  is  non-zero 
only  along  the  set  of  lines  defined  by 

w  =  ^a-t-Won  •  (2.6.66) 

Again,  the  set  of  offset  frequencies  may  be  obtained  by  evaluating  the  response  of  the 
system  to  a  constant  unit  input  sequence. 
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2.7  CHAPTER  SUMMARY 


In  this  chapter,  we  have  presented  some  of  the  fundamental  concepts  of  the  analysis 
of  linear,  discrete-time  time-varying  systems.  Both  time-domain  and  frequency-domain 
characterizations  of  discrete-time  LTV  systems  have  been  presented.  The  implications  of 
stability  and  causality  for  LTV  systems  have  been  discussed,  along  with  a  presentation  of 
the  analysis  of  LTV  network  structures.  The  discussion  of  classes  of  LTV  systems  will 
prove  valuable  to  us  when  we  begin  to  consider  modeling  discrete-time  nonstationary 
random  processes  using  LTV  systems. 


3.  MODELING  AND  FREQUENCY-DOMAIN  ANALYSIS  OF 
NONSTATIONARY  RANDOM  PROCESSES  VIA 
LINEAR  TIME-VARYING  SYSTEMS 


The  use  of  linear  time-invariant  system  models  has  proven  very  useful  in  the 
analysis  of  stationary  random  processes.  In  general,  one  may  assume  an  observed  wide- 
sense  stationary  random  process  is  the  output  of  a  LTI  system  excited  by  a  purely  random 
or  white  noise  process  (see  Koopmans,  1974,  pp.  233-238).  Frequency-domain 
characterizations  of  the  underlying  system  model  yield  insight  into  the  spectral  content  of 
the  observed  random  process  (see  Priestley,  1981,  pp.  267-268).  Time-domain 
characterizations  for  the  underlying  system  model  can  be  used  to  predict  values  of  the 
observed  random  process  or  to  filter  the  random  process  in  some  optimum  sense  (see 
Koopmans,  1974,  np.  212-217). 

In  a  similar  fashion,  linear  time-varying  system  models  are  useful  in  the  analysis  of 
nonstationary  random  processes.  In  this  chapter,  we  will  discuss  modeling  a  nonstationary 
random  process  as  the  output  of  a  LTV  system  excited  by  a  white  noise  process.  Spectral 
representations  of  the  nonstationary  random  process  which  can  be  derived  from  the 
frequency-domain  characterization  of  the  underlying  system  model  will  then  be  presented. 
A  discussion  of  the  implications  of  the  properties  of  the  observed  random  process  on  the 
characteristics  of  the  underlying  LTV  system  model  is  presented,  with  particular  attention 
paid  to  almost  cyclostationary  random  processes.  Finally,  the  characterization  of  cenain 
classes  of  nonstationary  random  processes  via  the  spectral  correlation  function  is  presented. 
For  a  more  thorough  treatment  of  these  topics,  see  Priestly,  1981,  pp.  817-866,  and 
Gardner,  1986,  pp.  301-359. 

3.1  LINEAR  TIME-VARYING  SYSTEM  MODELS  FOR 

NONSTATIONARY  RANDOM  PROCESSES 

In  order  to  introduce  the  concept  of  a  LTV  system  model  for  nonstationary  random 
proces.ses,  we  first  review  the  notion  of  a  LTI  system  model  for  stationary  random 
processes.  Let  {y(t)}  be  a  wide-sense  stationary  di.screte-time  random  process  such  that 

Ely{t)l  =  ()  (3.1.1) 


and 
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E[  1  y(t)  1^  ]  =  <  oo  . 

(3.1.2) 

By  the  Wold  decomposition  theorem  (see  Koopmans,  1974,  pp.  255-256),  we  know  that 
{y(t)}  can  be  expressed  as 

y(t) 

=  p(t)  -h  v(t) 

(3.1.3) 

where  the  processes  {|J.(t)}  and  {v(t)}  are  uncorrelated  zero-mean  wide-sense  stationary 
discrete-time  random  processes  with  the  following  properties: 

(i) 

{|i(t)}  can  be  represented  as  the  output  of  a  LTI  filter  excited 
by  a  wide-sense  stationary,  zero-mean  purely  random  (i.e., 
white  noise)  process.  In  other  words. 

oo 

^t(t)  =  X  h(i-X)  e(X) 

X  =  —  oo 

(3.1.4) 

where 

E[e(Xi)e*(X2)]=a^C(^i-X2)  . 

(3.1.5) 

Furthermore,  the  LTI  system  represented  by  h(X)  is  causal, 
and  has  the  property 

oo 

1  h(A.)  1  <  oo  . 

X  =  —  oo 

(3.1.6) 

(ii) 

The  process  (v(t)}  is  purely  deterministic,  in  other  words 
the  present  value  of  {v(t))  may  be  completely  determined  by 
linear  functions  of  its  past  values. 

For  the  moment,  we  will  consider  the  case  when  the  detemiinistic  component  of  {y(t)}  is 
not  present,  in  other  words  assume  v(t)  =  0  for  each  t.  We  make  this  assumption  without 
loss  of  generality  since,  from  property  (ii)  above  we  could  make  a  prediction  for  n(t)  from 
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past  observations  of  the  process  {y(t)}  and  remove  this  deterministic  component  from 
{y(t)).  We  may  therefore  without  loss  of  generality  express  {y(t)}  as 

oo 

y(0=  S  h(t-X)e(X)  .  (3.1.7) 

X  =  -  oo 

An  extension  of  the  Wold  decomposition  theory  for  nonstationary  discrete-time 
random  processes  is  given  in  Cramer,  1961.  Again  assume  that  {y(t)}  is  a  discrete-time 
random  process  satisfying 

E[y(t)]=0  (3.1.8) 

and 

E[ly(t)P]<°o  (3.1.9) 

for  each  t,  but  now  assume  {y(t))  is  nonstationary.  We  can  then  express  y(t)  as 

y(t)  =  p(t)  +  v(t)  (3.1.10) 

where  the  processes  {^.(t)}  and  {v(t)}  are  uncorrelated  zero-mean  discrete-time  random 
processes  with  the  following  properties: 

(i)  {lt(t)}  can  be  represented  as  the  output  of  a  LTV  filter 

excited  by  a  wide-sense  stationary,  zero-mean  purely 
random  process.  In  other  words, 

oo 

it(t)=  ]Lk(t,x)e(x)  (3.1.11) 

X  =  -  oo 

where 

E[t’(Xi)c*(X2)l  =  o‘'C(Xi-X2)  .  (3.1.12) 
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Furthermore,  the  LTV  system  represented  by  k(t,x)  is 
causal,  and  has  the  property 

oo 

X  1  k(t,t)  1^  <oo  .  (3.1.13) 

1=  -oo 

(ii)  The  process  {v(t))  is  purely  deterministic;  in  other  words, 
the  present  value  of  {v(t)}  may  be  completely  determined  by 
linear  functions  of  its  past  values. 

Consider  that,  since  v(t)  is  purely  deterministic,  we  can  predict  its  current  value  from  past 
observations  of  y(t)  and  so  remove  its  contribution  from  the  observed  random  process.  We 
may  therefore  without  loss  of  generality  restrict  our  attention  to  nonstationary  random 
processes  of  the  form 

oo 

y(t)=  Ek(t,t)e(T).  (3.1.14) 

X  =  -  oo 

Note  that  equation  (3.1.13)  does  not  ensure  BIBO  stability  for  the  LTV  system  for 
each  t.  However,  from  the  theory  of  Fourier  series  (see  Koopmans,  1974,  pp.  19-21)  we 
know  that  equation  (3.1.13)  holds  if  there  exists  some  function  K(t,a),  periodic  in  a  with 
period  2n,  which  is  continuous  for  almost  all  a  such  that 

oo 

K(t,a)  =  X  k(t,'c)  e**^  .  (3.1.15) 

X  =  -  oo 

Furthermore,  we  have  the  relationship 

7C  oo 

^  J  I  K(t,a)  1^  da  =  X  I  k(t,x)  1^  (3.1.16) 

-K  -oo 

for  each  t  (see  Koopmans,  1974,  p.  20).  We  will  make  use  of  the  fact  that  the  frequency 
response  function  for  the  LTV  system  model  exists  when  we  consider  frequency-domain 
repre.sentations  of  nonstationary  random  processes. 
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Notice  that  equation  (3.1.13)  ensures  that  the  LTV  system  model  will  exhibit 
limited  memory,  as  discussed  in  Section  2.2.  In  other  words,  the  contribution  to  some 
present  value  of  {y(t))  due  to  an  input  in  the  remote  past  is  negligible.  Causality  of  the 
LTV  system  model  guarantees  that  {y(t)}  does  not  depend  on  future  values  of  the  input. 
We  may  therefore  argue  that  input  values  contributing  significantly  to  {y(t)}  at  some 
present  time  are  distinct  from  those  input  values  contributing  significantly  to  {y(t))  in  the 
remote  past.  If  we  further  assume  that  the  {e(x))  are  independent  and  identically 
distributed,  then  v/c  may  make  a  mixing  assumption  for  {y(t)},  i.e.,  present  values  of 
{y(t)}  are  virtually  independent  of  values  of  {y(t))  in  the  remote  past.  This  mixing 
assumption  for  {y(t)}  will  ensure  that  appropriate  statistical  averages  involving  {y(t)}  are 
asymptotically  convergent  (see  Billingsley,  1968,  pp.  166-167). 

3.2  SPECTRAL  REPRESENTATIONS  FOR  NONSTATIONARY  RANDOM 
PROCESSES 

Assume  {y(t)}  is  a  nonstationary,  discrete-time  random  process  satisfying 
E(y(t)]=0  (3.2.1) 


and 


E[  ly(t)l^  ]  <  ~  .  (3.2.2) 

By  the  discussion  presented  in  the  previous  section,  we  know  that  y(t)  can  be  represented 
as 


oo 

y(t)-  S  k(t,T)e(t) 

T  =  -  oo 


(3.2.3) 


where  e(T)  is  a  zero-mean,  wide-.sense  stationary,  purely  random  discrete-time  process  with 
2 

variance  O  .  From  the  theory  of  spectral  repre.sentations  for  wide-sense  stationary  random 
prtK'esses.  we  know  that  (e(t))  can  be  represented  as 


E(t) 


lax 

e 


dZ(a)  . 


(3.2.4) 
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where  {Z(a))  is  a  zero-mean  complex- valued  random  process  over  (-7i,7c)  with  orthogonal 
increments;  in  other  words. 


E[dZ(ai)dZ*(a2)]  =  2no^  5(ai-a2)  dai  da2 


(3.2.5) 


(see  Priestley,  1981,  pp.  250-251). 


Now  consider  substituting  the  spectral  representation  for  e(T)  given  by 
equation  (3.2.4)  into  the  LTV  system  model  for  y(t)  given  by  equation  (3.2.3).  Our 
expression  for  y(t)  then  becomes 


oo  /  •• 

y(t)=  X  kit,-!)  ^  J 

X  =  -  oo  -7t 


lax  „ 

e  dZ(a) 


(3.2.6) 


Interchanging  the  order  of  summation  and  integration  in  equation  (3.2.6)  above,  we  obtain 


y(t)  = 


2k 


-7t 


/  oo 

1  I 


k(t.T)  e 


lat 


dZ(a) 


(3.2.7) 


Note  that  the  summation  in  brackets  in  (3.2.7)  above  is  the  frequency  response  function  for 
the  LTV  system,  K(t,a).  We  know  that  the  frequency  response  function  exists  from 
equation  (3.1.15).  Substituting  the  frequency  response  function  into  equation  (3.2.7) 
above,  we  obtain 


K 

y(t)  =  ^  j  K(t,a)  dZ(a)  .  (3.2.8) 

‘It 

Equation  (3.2.8)  defines  the  evolutionary  spectral  repre.sentation  of  the  nonstationary  (or 
evolutionary)  random  process  {y(t))  (see  Priestley,  1965,  and  Priestley,  1981,  pp.  825- 
826).  As  such,  it  is  analogous  to  the  spectral  representation  for  wide-sense  stationary 
random  processes.  We  will  see  that  this  representation  of  a  nonstationary  random  process 
is  useful  in  the  characterization  of  nonstationary  signals.  We  will  also  present  the 
implications  of  various  properties  of  the  prwess  (y(t) )  on  the  behavior  of  the  underlying 
system  model,  k(t,T). 
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3.3  CLASSES  OF  NONSTATIONARY  RANDOM  PROCESSES  AND 
IMPLICATIONS  OF  LINEAR  TIME-VARYING  SYSTEM  MODELS 


We  have  seen  in  the  previous  section  that  a  zero-mean,  finite  variance  nonstationary 
random  process  {y(t))  may  be  modeled  as  the  output  of  a  LTV  system  excited  by  a  wide- 
sense  stationary  random  process;  in  other  words, 

oo 

y(t)=  Xk(t,x)e(x)  .  (3.3.1) 

X  =  -  M 

The  evolutionary  spectral  representation  of  the  process,  given  by 
n 

y(i)  =  ^  j  K(t,a)  dZ(a)  ,  (3.3.2) 

‘7C 

was  also  presented.  Consider  from  equation  (3.3.2)  that  the  statistical  properties  of  the 
random  function,  {y(t)),  are  intimately  related  to  the  properties  of  the  frequency  response 
function  of  the  underlying  LTV  system  model,  K(t,a).  In  this  section  we  will  consider  the 
implications  of  the  system  response  function,  K(t,a),  on  the  statistical  properties  of  {y(t)} 
for  the  class  of  almost  cyclostationary  random  processes. 

Define  the  (symmetric)  autocorrelation  function  for  the  random  process  {y(t)}  as 

Ryy(a)  =  E  [  y(t+X)y*(t-)i)  ]  .  (3.3.3) 

The  random  process  {y(t))  is  said  to  be  almost  cyclostationary  if  Ryy(t,^)  is  an  almost 
periodic  function  of  t  for  each  X  (see  Gardner,  1986,  p.  302).  In  other  words, 

oo 

Ryy(t,M=  X  rn().)e“^"'  (3.3.4) 

n=:  oo 

where 
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(3.3.5) 


T 

,  lim  1  X-  ,  -id)„t 

rn(^)  =  XY  ^  Ryy(t.^)  e  ^ 

and  {())„}  is  a  sequence  of  increasing  frequencies  independent  of  X. 

The  random  process  { y(t) }  is  almost  cyclostationary  if  the  underlying  LTV  model 
for  {y(t))  has  a  frequency  response  function  which  is  almost  periodic  in  t  for  all  a,  in  other 
words  if  k(t,T)  represents  a  LPTV  system.  In  order  to  understand  this,  first  assume  that  the 
frequency  response  function  for  the  underlying  LTV  system  model  is  almost  periodic  in  t. 
From  the  discussion  of  almost  periodic  frequency  response  functions  presented  in  Chapter 
2,  we  know  that 


K(t,a) 


Z  an(a) 


io)n(a)t 

e 


n=  -«> 


(3.3.6) 


where 


lim 


k) 


=  ^  K:(t,a)  e . 

to-^OO  2to 


(3.3.7) 


From  the  LTV  system  model  for  {y(t))  given  by  (3.3.1),  we  obtain 


oo  71  ‘ ♦ 

yd)  =  S  ^  I  an(a)  e‘  "  dZ(a)  . 

n=  -oo  -K 


(3.3.8) 


Substituting  equation  (3.3.8)  into  the  expression  for  Ryy(t,^)  given  by  equation  (3.3.3) 
and  simplifying,  we  obtain 


Ryy(t,>.)  =  X  k(t+>„x)k*(t-A.,x)  . 


(3.3.9) 


X=  -oo 


It  is  straightforward  to  show  that 
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2  n 

Ryy(t,X)  =  ^  J  )c(t+X.a)K:*(t-X.a)  da  . 


(3.3.10) 


Substituting  for  K(t,a)  from  equation  (3.3.6)  into  (3.3.10),  we  obtain 


oo  oo 


Ryy(t,A.)  =  ^  J  S  S  an(a)a*m(a)  e 
-Jt  n=  m=  — oo 


*  ia)„(a)(t+A.)  -ico^(a)(t-X) 

m(a)  e  "  e  da . 


(3.3.11) 

Recall  that  the  LTV  system  represented  by  k(t,t)  is  causal,  and  so  from  the  discussion  of 
causal  LPTV  systems  presented  in  Chapter  2,  we  know  that 


a)n(a)  =  pa  +  tOo 


(3.3.12) 


Our  expression  for  Ryy(t,X)  can  therefore  be  rewritten  as 


2  oo  oo 


Ryy(t,A.)  =  2^  ^  2^  J  3fi(0i)j  in(®)®  ® 


n=— oo  ni=-oo  -K 


Notice  that  the  sequence 


(3.3.13) 


{O)on-Wom:  n  =  -I.  0.  1.  •••  ;  m  =  ....  -1,  0,  1, 


(3.3.14) 


can  be  thought  of  as  a  sequence 


|(}>j:  j  =  ..., -1,0,  1,...} 


(3.3.15) 


of  distinct,  increasing  frequencies.  Now  define 


PjlA)  —  _  y ^  Jjp(oc)jpi(oc;L  c 


NMj  -n 


(3.3.16) 


where  the  set 
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NMj  —  {  n,m.  COon't*Joni  ~  ) 

(3.3.17) 

Our  expression  for  Ryy(t,X,)  may  therefore  be  rewritten  as 

oo 

Ryy(t,X)  =  X  r,a)  ; 

j=-oo 

(3.3.18) 

in  other  words,  Ryy(t,A.)  is  an  almost  periodic  function  in  t.  So,  if  the  underlying  LTV 

system  frequency  response  function  is  an  almost  periodic  function  of  t,  then  {y(t)}  will  be  ‘ 

an  almost  cyclostationary  process. 

Next,  assume  that  {y(t)}  is  almost  cyclostationary,  but  the  frequency  response 
function  for  the  underlying  LTV  system  model  is  not  almost  periodic.  In  other  words, 

assume  that 

lirn  1  -itot 

,  ^  ^  L  K(t.a)  e  =0 

>oo  2to 

(3.3.19) 

for  any  choice  of  to  or  a.  Consider  that  equation  (3.3. 19)  implies  that 

oo 

-itot 

2^  K(t,a)  e 

t=-oo 

(3.3.20) 

exists  for  all  O)  and  a.  Since  the  product  of  convergent  series 
equation  (3.3.20)  implies  that 

is  convergent. 

oo  oo 

V*  V'  /  V  -iOiti  -ici)2t2 

X  2^  K(ti,a)K:  (t2,a)  e  ‘‘e 

t]  =  -oo  l2  = 

(3.3.21) 

exists,  or  by  a  change  of  vanables. 

-IWA. 

X  X  ^■U+/.,a)K'  (t-/.,a)  e  e 

A.  =  -oo  1  =  -oo 

(3.3.22) 
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exists  for  all  <|>  and  Vj/. 


Since  {y(t)}  is  almost  cyclostationary. 


lim  1  ^ 

T-><»2T  Ryy(t.A.) 


(3.3.23) 


exists  and  is  nonzero  for  some  sequence  of  {(j),,)  and  some  X.  From  equation  (3.3.10),  we 
obtain 


n 


r 

J 


-K 


lim 

T->o. 


1 

2T 


X  K(t+>.,a)K*(t-X,a)  e 
i=-T 


da 


(3.3.24) 


exists  and  is  nonzero.  Notice  that  equation  (3.3.24)  implies  that 


lim  1 
T->oo  ZT 


^  X  K:(t+>.,a)K*(t-X,a)  e 


(3.3.25) 


t=-T 


exists  and  is  nonzero  over  some  interval  in  a.  Moreover,  equation  (3.3.25)  implies  that 


^  K(t+X,a)K*(t-X,a)  e 


(3.3.26) 


t=-oo 


diverges.  Notice  that  if  the  sum  in  equation  (3.3.26)  diverges,  then  the  sum  in 
equation  (3.3.22)  also  diverges,  which  is  a  contradiction. 

From  the  discussion  presented  above,  we  see  that  the  characterization  of  almost 
eyelostationary  random  processes  is  intimately  related  to  the  characterization  of  LTV 
systems  with  almost  periodic  frequency  response  functions.  We  will  make  use  of  this  fact 
when  we  consider  the  empirical  analysis  of  LTV  systems  in  Chapter  4, 
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3.4  SPECTRAL  CORRELATION  FUNCTIONS  FOR  NONSTATIONARY 
RANDOM  PROCESSES 

The  conventional  power  spectrum  has  proven  to  be  a  useful  means  of  characterizing 
wide-sense  stationary  random  processes  in  the  frequency  domain.  Varicus  extensions  of 
the  concept  of  the  power  spectrum  for  nonstationary  random  p;.jcesses  have  been  proposed 
over  the  years  (see  Lyones,  1968,  and  Cohen,  1989,  for  a  review  of  the  subject).  In  this 
section,  we  will  discuss  a  technique  which  has  proven  useful  in  the  study  of  finite-energy 
and  almost  cyclostationary  random  processes.  This  technique  is  based  on  the  fact  that,  for 
a  nonstationary  random  process,  correlation  exists  between  distinct  frequency  components 
of  the  spectral  representation.  Estimates  of  the  spectral  correlation  function  therefore 
provide  a  means  of  characterizing  the  random  process  in  the  frequency  domain.  For  a  more 
thorough  treatment  of  the  subject  of  spectral  correlation  functions,  see  Gardner,  1986, 
pp.  .^01-359,  and  Bendat  and  Piersol,  1986,  pp.  447-468. 

Assume  that  {x('C))  and  {y(t)}  are  discrete-time  random  processes  such  that 


E(x(t)1=0,  (3.4.1) 

E(  I  x(T)  r  1  <oo  ,  (3.4.2) 

and 

lx(x)l<c«  (3.4.3) 

for  each  x  and  for  every  realization  of  (x(t)),  and  such  that 

E[y(t)|=().  (3.4.4) 

El  I  y(t)  r  1  <«,  ,  (3.4.5) 


and 

ly(t)l<c«  (3.4.6) 
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for  each  t  and  for  every  realization  of  {y(t)}.  Define  the  finite-duration  Fourier  transforms 
of  {x(t)}  and  {y(t)}  as 


to 

Xo(a)=  X  x(T)e'‘“^  (3.4.7) 

t=-to 


and 


k) 

Yo(0))  =  X  y(0  e  .  (3.4.8) 

t=-to 

Note  that  both  Xo(a)  and  Yo((o)  exist  since,  from  equations  (3.4.3)  and  (3.4.6),  they  are 
finite  sums  of  finite-amplitude  sequences.  Formally  define  the  cross-spectral  correlation 
function  as 

SYx(to,a)  =  E  [Yo(a))Xo*(a)]  (3.4.9) 

(see  B'.ndat  and  Piersol,  1986,  pp.  448-450).  Notice  that,  in  some  sense,  the  cross- 
spectral  correlation  function  indicates  the  amount  of  correlation  between  distinct  spectral 
components  of  the  two  random  processes,  {x(t))  and  {y(t)).  Substituting  the  definitions 
for  the  finite-duration  Fourier  transforms  given  by  equations  (3.4.7)  and  (3.4.8)  into 
(3.4.9),  we  obtain 


'to 

Ziax  -1(1)1 

RYx(t.x)e  e  (3.4.10) 

i:=-To 

where  we  have  defined  the  (asymmetric)  autocorrelation  function 

RYx(t,T)  =  E[y(t)  x*(x)l  .  (3.4.11) 

Assuming  that  the  limit  in  equation  (3.4.10)  exists,  we  obtain 


SYx(0),a)  = 


lim  lim 

to— )cx3  tq— 


to 

I 

i=-to 
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SYx(t>J'Ct)  =  X 

t=-oo 


I 


RYx(t.'^) 


iax  -itot 
e  e 


Xs-oo 


(3.4.12) 


For  the  moment,  assume  that  {x(t))  and  {y(t)}  are  wide-sense  stationary  random 
processes  so  that 

RYx(t.t)  =  %x(t-t)  .  (3.4.13) 


From  equation  (3.4.12),  we  obtain 


icxx  -icoi 

SYx(w,a)  =2-  2.  %x(t-t)  e  e 


(3.4.14) 


or  by  the  change  of  variables  X  =  t-x. 


SYx(w.a) 


X=-oo 


X  X  ^Yx(^)  e 


i(oX 


h=- 


-i(o)-a)x 


(3.4.15) 


Notice  that  the  summation  in  brackets  in  equation  (3.4.15)  above  is  the  conventional  cross¬ 
power  spectrum,  Syx(^)  (see  Bendat  and  Piersol,  1986,  p.  121).  Equation  (3.4.15) 

therefore  becomes 

SYx(0J,a)  =  SYx(tt))  5(o>-a)  .  (3.4.16) 


In  other  words,  wide-sen.se  stationary  random  processes  exhibit  spectral  correlation  only 
on  the  line  (o  =  a  in  the  (co,a)  bifrequency  plane. 


Consider  the  system  indicated  schematically  in  Figure  3.4.1.  The  purely  random 
process  {e(x)}  excites  an  LTV  system  represented  by  k(t,x)  to  produce  the  nonstationary 
random  process  {y(t)}.  This  .same  purely  random  process  also  excites  the  LTI  system 
represented  by  h(x)  to  produce  the  wide-sense  stationary  random  process  {x(x)).  From  the 
discussion  of  Chapter  2,  we  know  that  x(x)  is  given  by 
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OP 

X(x)=  X  h(x-X)e(X)  .  (3.4.17) 

X=.oo 


From  the  spectral  representation  of  {e(x))  given  by  equation  (3.2.4),  we  obtain 


X  h(x-X)  e‘^^  dZ(7)  . 


or  simply 


71 


H(y)  e‘^  dZ(y)  , 


(3.4.18) 


(3.4.19) 


where  H(y)  is  the  frequency  response  function  of  the  LTl  system  given  by 
equation  (2.3.6).  Substituting  equation  (3.4.19)  into  the  definition  of  Xo(a)  given  by 
equation  (3.4.7),  we  obtain 


Xo(a)  =  ^  J 


Vx=-xo 


(3.4.20) 


Recall  from  equation  (3.2.8)  that  {y(t))  admits  the  evolutionary  spectral  representation 


n 


y(t)  = 


J_ 

2rc 
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K(t,a)  dZ(a)  . 


(3.4.21) 


Substituting  equation  (3.4.21)  into  the  definition  for  Yo(w)  given  by  equation  (3.4.8),  we 
obtain 


K 

-7t 


z 


V‘=-'o 


K(t,a) 


-itoi 

e 


IdZ(a)  . 


J 


(3.4.22) 
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From  the  definition  of  the  cross-spectral  correlation  function  given  by  equation  (3.4.9),  we 
obtain 


lim  lim  r  *  t 

Svx((i).a) 

7C 

=  ^  J  K(a),p)  H*({3)  6(i3-a)  dP  ,  (3.4.23) 

-K 

or,  simply, 

SYx(0),a)  =  ^  K(co,a)  H*(a)  .  (3.4.24) 

Notice  that  equation  (3.4.24)  gives  a  very  straightforward  relationship  between  the  cross- 
spectral  correlation  function  for  {y(t)}  and  {x(t)}  and  the  linear  systems  which  produced 
these  random  processes.  We  will  make  use  of  this  fact  when  we  consider  the  empirical 
analysis  of  linear  time-varying  systems  in  Chapter  4. 


AS-9 1-367 


FIGURE  3.4.1 

LTV  AND  LTI  SYSTEMS  EXCITED  BY  A  PURELY  RANDOM  PROCESS 
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3.5  CHAPTER  SUMMARY 


In  this  chapter  we  have  introduced  some  of  the  fundamental  concepts  regarding  the 
modeling  of  nonstationary  random  processes  using  linear  time-varying  filters.  The  use  of  a 
frequency-domain  characterization  of  the  LTV  system  model  to  obtain  an  evolutionary 
spectral  representation  for  a  nonstationary  random  process  was  also  presented.  Finally,  the 
relationship  between  the  LTV  system  models  for  nonstationary  random  processes  and  the 
spectral  correlation  function  for  these  random  processes  was  presented.  In  the  following 
chapters,  these  concepts  will  be  applied  to  the  characterization  of  LTV  systems  excited  by 
random  signals. 


.*>4 


4.  CHARACTERIZATION  OF  LINEAR  PERIODICALLY 
TIME-VARYING  SYSTEMS 


In  previous  chapters,  we  reviewed  the  principal  methods  for  analytically  describing 
a  discrete  time  linear  time-varying  system  and  noted  the  properties  of  such  system 
descriptions.  Techniques  for  the  frequency-domain  analysis  of  the  nonstationary  signals 
arising  as  the  outputs  of  time-varying  systems  were  also  presented.  We  now  turn  our 
attention  to  the  problem  of  empirically  characterizing  LPTV  systems.  To  be  more  specific, 
we  require  techniques  for  obtaining  a  reliable  estimate  of  the  system  transmission  function, 
K(o),a),  for  a  given  LPTV  system.  The  development  and  assessment  of  these  techniques 
will  rely  heavily  on  the  concepts  and  definitions  presented  earlier. 

4.1  CHARACTERIZATION  OF  LINEAR  PERIODICALLY  TIME- 

VARYING  SYSTEMS  via  DETERMINISTIC  SIGNALS 

Recall  the  definition  of  the  time-varying  system  transmission  function,  K(a),a), 
given  by  equation  (2.3.25) 

oo 

K((o,a)  =  X  K(t,a)  e ,  (4.1.1) 

I  =  -  oo 

where  K(t,a)  is  the  time-varying  frequency  response  function.  Since  K(t,a)  is  defined  as 
the  response  of  the  LPTV  system  to  a  complex  sinusoid  of  frequency  a,  the  system 
transmission  function  can  be  viewed  as  the  Fourier  transform  of  the  output  signal  produced 
by  an  applied  complex  sinusoid  of  frequency  a.  This  observation  suggests  that  one 
straightforward  method  of  estimating  the  system  transmission  function  for  a  given  LPTV 
system  would  be  to  apply  a  complex  sinusoid  of  a  known  frequency,  observe  the  resulting 
output  signal,  and  compute  an  estimate  of  its  Fourier  transform.  This  would  yield  an 
estimate  of  the  system  transmission  function,  K(co.a),  for  some  specific  input  frequency, 
a.  Estimates  of  the  system  transmission  function  for  other  input  frequencies  could  be 
obtained  by  repeating  this  process  for  each  desired  input  frequency. 

Consider  applying  a  complex  sinusoid  of  freciuency  a  to  a  causal,  stable  LPTV 
system.  Assume  that  the  input  signal,  x(T),  as  given  by 

KIT 


\(x)  =  costal)  +  i  sintax)  =  e 


(4.1.2) 


has  been  applied  to  the  system  of  interest  for  a  sufficient  time  for  the  system  to  have 
reached  a  steady-state  condition.  If  we  obtain  N  samples  of  the  resulting  output  sequence, 
y(t),  we  may  compute  the  Fourier  transform  of  the  output  sequence,  Y(to),  as 

N-l 

Y((o)=  X  K(t,a)  e .  (4.1.3) 

t  =  0 


From  the  propenies  of  the  Fourier  transform,  we  obtain 
n 

Y(to)  =  ^  I  K(Y,a)  FnCy-O))  dy  , 


where  Fn(y)  is  the  Fejer  kernel,  defined  as 
N-l 

Fn(y)=  I  e-’" 

t  =  0 


(4.1.4) 


(4.1.5) 


(see  Papoulis,  1962,  pp.  42-47).  For  N  sufficiently  large,  the  Fejer  kernel  may  be 
approximated  by  a  Dirac  delta  function,  and  so 

Y(a))  =  K(o),a)  .  (4.1.6) 

Several  points  can  be  made  concerning  the  method  of  estimating  K(a),a)  presented 
above.  First,  the  accuracy  of  the  approximation  of  K(a),a)  by  Y(a))  is  intimately  related  to 
the  available  length  of  observed  data,  N.  This  issue  of  distortion  of  the  estimate  of  K((D,a) 
due  to  finite  data  lengths  closely  parallels  the  issue  of  distortion  of  spectral  estimates  due  to 
finite  data  lengths.  The  familiar  method  of  applying  a  windowing  function  to  the  observed 
data  to  reduce  the  effects  of  finite  data  lengths  could  also  be  applied  to  improve  our  estimate 
of  K(a),a)  (see  Harris,  I97S).  Secondly,  this  method  provides  an  estimate  for  K(o),a) 
for  some  specific  input  frequency,  a.  Estimation  of  K(to,a)  for  other  values  of  a  requires 
that  the  frequency  of  the  applied  signal  be  changed  and  that  a  new  record  of  the  output 
sequence  be  obtained  after  sufficient  time  has  passed  for  a  steady-state  condition  to  be 
reached. 
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The  final  point  to  be  made  regarding  this  method  relates  to  the  generation  of  the 
applied  input  signal.  For  a  digital  implementation  of  the  LPTV  system,  obtaining  a 
complex  sinusoidal  input  is  a  relatively  straightforward  matter,  involving  no  more  than  the 
generation  of  a  sequence  of  sines  and  cosines.  However,  in  a  LPTV  system  in  which 
analog  signals  are  digitized  and  then  processed  to  obtain  digital  sequences,  it  is  impossible 
to  generate  a  complex  input  signal.  One  may  overcome  this  problem  by  applying  two 
separate  input  signals,  cos(aT)  and  sin(at),  and  computing  the  Fourier  transforms  for  each 
resulting  output  sequence  separately.  An  estimate  for  K(a),a)  could  then  be  obtained  from 
these  two  Fourier  transforms.  To  illustrate,  let  yR(t)  and  yi(t)  designate  the  outputs 
obtained  from  the  LPTV  system  due  to  application  of  the  input  signals  cos(aT)  and 
sin(ax),  respectively.  By  linearity  of  the  system  of  interest  and  of  the  Fourier  transform, 
we  have 

YrIco)  +  i Y,((o)  =  Y((o)  =  K((o,a)  ,  (4.1 .7) 

where  YR(a))  and  Yi(to)  are  the  N-length  Fourier  transforms  of  yR(t)  and  yi(t), 
respectively. 

Note  that  the  accuracy  of  our  estimate  for  the  system  transmission  function  is  now 
dependent  on  both  the  length  of  the  observation  interval  and  on  the  relative  phase  of  the  two 
input  signals.  Assume  that  there  is  a  phase  error  of  ())  between  the  two  input  signals;  in 
other  words,  assume  that  we  apply  cos(ax)  and  sin(ax  +  0)  at  the  input  of  the  LPTV 
system.  For  small  values  of  phase  error, 

sin(ax+(j))  =  sin(ax)  +  <|)cos(ax)  ,  (4.1.8) 


so  that 


xo(x)  =  cos(ax)  +  1  sin(ax+0)  =  e  +9^  (4.1.9) 

l-omiing  our  estimate  for  the  system  transmi.ssion  function  as  before,  we  obtain 

Y0(o))  =  K(co,a)  +  ^  K((o.-a)  .  (4.1.10) 
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The  significance  of  this  distortion  in  our  estimate  of  the  system  transmission  function 
depends  on  the  magnitude  of  the  phase  error  and  on  the  nature  of  K(a),a). 

4.2  CHARACTERIZATION  OF  LINEAR  TIME- VARYING  SYSTEMS  via 
NONDETERMINISTIC  SIGNALS 

We  now  consider  a  technique  for  estimating  the  system  transmission  function  of  a 
LKrV  system  by  observing  the  response  of  the  system  to  a  nondeterministic  input  signal. 
This  technique  is  based  on  the  observation  made  in  Section  3.4  that  the  cross-spectral 
correlation  function  between  the  output  of  a  LPTV  system  and  the  random  input  to  the 
system  is  directly  proportional  to  the  LPTV  system  transmission  function.  In  this  section, 
we  will  derive  the  approximate  statistical  properties  of  the  estimated  cross-spectral 
correlation  function  and  of  the  estimated  .system  transmission  function. 

Let  {e(x))  by  a  wide-sense  stationary,  complex  Gaussian  purely  random  process 
satisfying 

E(e(T)l=0  ,  (4.2.1) 


and 


E  (  e(ti)e*(X2) 


o2  ;  Xi  =  X2 
0  ;  else 


(4.2.2) 


Next,  assume  that  (e(x)}  is  applied  to  the  input  of  a  causal,  stable  LPTV  system  with  the 
simple  .system  transmission  function 

K(a),a)  =  Ho(a)  6(o)-(o„  -  ~  >  (4.2.3) 

where  p  is  the  rate  of  the  system.  Note  that  for  this  LPTV  system,  K(a),a)  is  non-zero 
only  along  the  line 


to  =  co^, 


(4.2.4) 


More  sophisticated  LFPV  systems  may  be  considered  to  be  linear  combinations  of  systems 
of  this  ptirticular  fomi. 


Recall  from  equation  (3.2.4)  that  the  spectral  representation  for  {£(x))  is  given  by 


Ik 


e^  dZ(P) 


(4.2.5) 


The  K-length  discrete  Fourier  transform  of  (e(x)}  is  given  by 


£K(ak)  =  Z  e('t)  e  . 

T  =  0 


(4.2.6) 


where 


27tk 

ak  =  ^,  k  =  0,  1....  .K-1 


(4.2.7) 


Substituting  for  the  spectral  representation  for  {e(x)}  into  equation  (4.2.6),  we  obtain 


^Klotk)  =  ^  J  FK(P-ak)  dZ(p)  , 


(4.2.8) 


where  Fk(P)  is  the  Fejer  kernel  of  equation  (4.1.5).  From  equation  (4.2.8)  we  see  that  the 
{  £K(ctk) )  are  a  set  of  zero-mean,  complex  Gaussian  random  variables.  Furthermore, 


note  that 


1  II 

E  £  k(«j)  J  FK(P'ak)  F*k(P-cCj) 


(4.2.9) 


or  simplifying. 


E  £K(ttk)  £  K(OIj)  = 


i  K  a-  ;  k  =  j 
1  0  ;  else 


(4.2.10) 
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so  that  {  £K(ctk) }  are  a  set  of  i.i.d.  complex  Gaussian  random  variables  with  zero  mean 
and  variance  Ka^. 

Now  define  the  N-length  discrete  Fourier  transform  of  the  system  output,  y(t),  as 


N-i 


YN(O)n)  =  X  y(0  e 

1=0 


(4.2.11) 


Recall  from  equation  (3.2.8)  that  we  may  express  y(t)  as 


y(t)  =  ^  J  K:(t,a)dZ(a)  , 

-7t 


(4.2.12) 


where  K(t,a)  is  the  system  response  function.  For  the  LPTV  system  under  discussion. 


1  r 

K(t,Y)  =  ^  J  K(0),y)  e  dt  , 


(4.2.13) 


K(t,Y)  =  Ho(Y)  exp{  i(a)„ +^)t  )  .  (4.2.14) 

Substituting  for  the  system  response  function  in  equation  (4.2.12),  we  obtain 

K 

y(t)  =  ^  J  i(a)o+^)t  )  dZ(Y)  ,  (4.2.15) 

-7t  P 


so  that 


YN(Cl>n)  = 


r,  I 


Fn  (Wn-tOo  - 


1 

P 


)  dZ(Y) 


(4.2.16) 
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Notice  that  the  {YN(o)n))  are  a  set  of  complex,  zero-mean  Gaussian  random  variables. 
Furthermore,  consider  that 

jt 

E  { YN(a)n)  Y*n(o)j)}  j  lHo(Y)l^  Fn  ((0„-(0o -h  F*n  (cOj-tOo  -h  dY(4.2.17) 

-Jt  P  P 


so  that  for  N  sufficiently  large, 

E  {YN(to„)  yVo)j)1  =  I  lHo(p(a)„-(Oo))l^  n  =  j  (4.2.18) 

10  ;  else 

Now  consider  choosing  N,  the  length  of  the  output  DFT,  and  K,  the  length  of  the 
input  DFT,  so  that 

N 

K  =  P-  (4.2.19) 


c  27tn  .  27ik 

For  any  (o,,  =  and  any  we  may  write 


(On  = 


27i(n-k)  27tk 


N  pK  ~  p 


(4.2.20) 


From  equation  (4.2.20)  we  note  that  the  set  of  points 

{  (^l^n+tn'0(k+m)'  ~  i  M/2)  (4.2.21) 

lie  along  a  line  of  slope  -  passing  through  the  point  (cOn.ctk)- 


Now  consider  the  two-dimensional  complex  random  variable 


1  o* 

KtOn.ttk)  =  ^  YN((On)  C  KlOCk)  • 


(4.2.22) 


Note  that 
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(4.2.23) 


E  {l(0)n,ak)}  =  i  ^  I  lijCY)  FN(t«)„-ti)o-^)  F*K(Y-ak)  dY  . 

-71  ” 


or,  for  N  and  K  sufficiently  large, 


E  {I(cOn,ak)}  = 


a2Ho(p(C0n-0)o))  ;  Wn  =  tOo  + 


ttk 


(4.2.24) 


0 


;  else 


Furthermore,  note  that 

E  {  I  I(o)n,ak)  1^}  =  I  E  {l(C0n,ak)}  l^+yl  Ho(p(cOn-(Oo))  1^  ,  (4.2.25) 


so  that 


Var  {  l(cOn,ak)  )  =  ^  *  Ho(p(a)n-C0o))  1^  ■  (4.2.26) 

From  equations  (4.2.24)  and  (4.2.26)  we  see  that  I  (o>n,ak)  is  an  inconsistent  estimator 
for  the  system  transmission  function. 

Recall  that  our  choice  of  N  and  K  led  to  a  set  of  discrete  frequency  points  such  that 
{  ((On+m,ak+m):  HI  =  0,  ±1,  ...  ,  ±  M/2)  lies  along  a  line  of  slope  1/p  in  the  (0),a) 
bifrequency  plane.  If  we  assume  that  K(o),a)  is  essentially  constant  over  this  set  of  points, 
we  could  consider  estimating  K(co,a)  by  smoothing  I(t0n,ak)-  Define  the  random  variable 

M/2 

Sm  ("'nittk)  riTT  ^  ^  (Wn+m»ttk+m)  »  (4.2.27) 

m  =  -M/2 


and  note  that 


M/2 

E  {  SMtCOfi.Ctk)  }  =  I  ^  E  {  l((l)n+in'®k-'-m)  ) 

m=-M/2 


(4.2.28) 
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For  N  and  K  sufficiently  large,  and  assuming  that  K((0,a)  is  constant  over  the  set  of  points 
{  (0)n+m.ctk+m):  m  =  0,  ±1, ,  ±  M/2  ).  we  obtain 


E  {  Sf4(C0n»Cilt)  }  ““  1 


o^Hofttk) 

0 


Ctk 

;  0)n  =  too  +  y 
;  else 


(4.2.29) 


Further,  note  that 

E{  I  SM((0n,ak)  l'}  =  I  E  (SM(a)n.ak))  +  I  Ho(p(a)„-C0o))  1^  ,  (4.2.30) 


so  that 


Var  {  SM(Wn,ak)  }  =  •  Ho(p((0n-t0o))  1^  .  (4.2.31) 

p(M  +  i ) 

If  we  pick  M  so  that  as  M,  N,  and  K  increase  without  bounds  M/N  0  and  M/K  — >  0, 
then  SM(C0n,ak)  is  a  consistent  estimator  of  the  system  transmission  function. 

We  now  consider  the  statistical  properties  of  SmIWik^Ic)-  First,  notice  from 
equation  (4.2.27)  that 


M/2 


Sm  ((On,ak)  = 


M+i 


+m  )£  K(CCk  +m) 

m  =  -M/2 


(4.2.32) 


For  M  sufficiently  large,  and  for  tOn  ^  O),,  + 


ttk 


(  YN(Wn^m)£  K(ak+m);  m  =  0,  ±  1 ,  ,  ±  M/2  j 


(4.2.33) 


o 


4 


are  a  set  of  uncorrelated,  identically  distributed  complex  random  variables.  We  can 
therefore  argue  that  SM(()\,.ak)  is  approximately  Gaussian  distributed  with  zero  mean  and 

For  oon  =  o)„  -1-  — \  SMtcOn.ttk)  is  approximately 

2 

Gaussian  distributed  with  mean  o^H„(p(tOn-oj„))  and  variance  I  Ho(p(t0n‘Wo))  I  • 

p(M  +  l) 


var’ nice - i  H„( p(a)n-to„) ) 

p(M‘l) 
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Finally,  assume  that  we  estimate  the  noise  variance,  0^,  from  the  given 
observations  of  {e(x))  as 


x=0 


(4.2.34) 


where  %  2K  is  a  chi-square  distributed  random  variable  with  2K  degrees  of  freedom.  Our 
estimate  for  the  system  transmission  function  is  therefore  given  by 


^((On.CXk)  —  A  'SM(tOn,(X|c) 
^2 


(4.2.35) 


For  K  large,  02  is  approximated  by 


^2  ~  o2  n 


(4.2.36) 


where  t)  is  a  Gaussian  random  variable  with  unit  mean  and  variance  1/K.  Our  estimated 
system  transnrission  function  is  therefore  approximated  by 


K(C0n,ak)  = 


K(C0n,ak)  +  — 


fp(M+l)  ^ 


Ho(p(C0n-{0o)) 


fp(M+l) 


p  ;  else 


(4.2.37) 


where  p  is  a  complex,  zero-mean,  unit  variance  Gaussian  random  variable. 

Now  consider  a  more  sophisticated  LPTV  system  composed  of  linear  combinations 
of  systems  of  the  fomi  assumed  above.  In  other  words,  assume  that 


K{to,a)  =  x  Ho(a)  5((o-w„- -  ) 
o  P 


(4.2.38) 


for  some  set  of  distinct  (o)„l.  Notice  that  SkIocr)  ■''tiH  de.scribed  by  equation  (4.2.8), 


but  that  now 


n 

YN(a)n)  =  Z  Ho(Y)  FN((0„-a)o-^)  dZ(Y)  .  (4.2.39) 

o  -n  ” 

Since  the  (Wi,)  are  distinct,  for  N  sufficiently  large  we  obtain 

■  ^  2 

E  {  YN(a)n)  Y*n(C0j)  )  =  ^  o  ;  n=j  (4.2.40) 

-  0  ;  e 1 se 


If  we  again  define  I(tOn,ai()  by  equation  (4.2.23),  we  obtain 

E  {  KtOn.ttk)  }  -<  <y^Ho(ak)  ;  cOn  =  cOo  +  y  (4.2.41) 

0  ;  else 

and 

Var  {  I(o)n.ak)  )  =  Ho(p(a)„-a)o))  1^  .  (4.2.42) 

P  0 

Using  equation  (4.2.28)  as  our  definition  for  SM(tOn.ctk).  we  obtain 

E  {  SM(cOn,ak)  }  =  ‘  f^^Ho(ak)  ;  ti)„  =  tOo  +  -^  (4.2.43) 

0  ;  else 

and 


Var  {  SM(C0n,ak)  }  jX  H„(p(0)n-(0o))  1^  •  (4.2.44) 

Defining  our  estimate  for  the  system  transmission  function  by  equation  (4.2.35)  yields  a 
complex  random  variable  which  is  approximated  by 
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X  Ho(p(0)„-C0o)) 


i^(Wn,ak)  =  ^ 


K(C0n,ak)  + 


Vp(M+l) 

X  Ho(p(0)„-COo)) 
_o _ 

Vp(M+l) 


p  ;  Q)„  =  a)o+  ^ 


,  (4.2.45) 


p  ;  else 


where  p  is  a  complex  Gaussian  random  variable  with  zero  mean  and  unit  variance. 


Finally,  consider  the  case  of  applying  this  technique  for  estimating  K(a),a)  when 
both  {e(t))  and  {y(t)}  are  observed  in  the  presence  of  additive,  stationary  independent 
white  noise,  as  indicated  schematically  in  Figure  4.2.1.  Assume  that  the  input  observation 
noise,  {<t)(x)),  and  the  output  observation  noise,  {^(0).  are  wide  sense  stationary  complex 
Gaussian  random  variables  such  that 


E  { 


C(tl)C*(t2)  }  = 


j  0^2  ;  ti  =  t2 
[0  ;  else 


and 


E  {  <)(t,)(i>*(X2)  }  = 


Xl  =  T2 
else 


(4.2.46) 


(4.2.47) 


Further,  assume  that  {({)('t)}  and  {^(t)}  are  mutually  independent,  and  independent  of 
{e(T)}. 


The  discrete  Fourier  transforms  of  the  observed  input,  {r(x)),  and  the  observed 
output,  {x(t)),  are  given  by 


K-1  y, 

RKlctk)  =  X  •■(X)  e  =  £K(ak)  +  CK(ak). 

T  =  0 


(4.2.48) 


and 
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(4.2.49) 


N-1 

XN(0)n)=  2  X(t) 
t  =  0 


e  =  Yn(CO„)  +  ON(tOn) 


where 


_L 

2n 


J 


FK(Y-ak)  dZ;(Y) 


and 


N-i 

^N(o)n)  =  X  <t>(0  e 

t  =  0 


-‘“nl 


=  ^  I  FK(Y-ak)  dZ0(Y) 


I 


-7t 


Now  define  the  random  variable  I'CoJn.ak)  as 


r(0)n,ak)  =  ^  XN(COn)  R*K(ak)  . 


and  let 


M/2 

S  M(^n’®k)  ~  jvTiT  ^  (^n+m»®k+m)  » 

m=-M/2 

SO  that 

M/2  ^ 

S  M(t^n’®k)  M+1  ^  |S(  YN(tt>n+m)  ^  K(®k+m) 

m=-M/2 

M/2  ^ 

M  +  1  ^  N  ^N(^+m)  K(®k+in) 

m=-M/2 

1  1  r  * 

^  M+l  ^  N  ^N(^+m)  ^  K(Ctk+m) 
m=-M/2 

M+i  w  ^N(Wn+m)  K(Otk+m) 

m=-M/2 


(4.2.50) 


(4.2.51) 


(4.2.52) 


(4.2.53) 


(4.2.54) 
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Consider  the  term 


1 

M+l 


M/2  ^  * 

X  'n  ^N(t^+m)  £  K(Otk  +fn) 


in=-M/2 


(4.2.55) 


in  equation  (4.2.55)  above.  Since  {^(t)}  and  {e(t)}  are  independent  Gaussian  random 
variables. 


E  {  N  ^  K(ak)  }  =  0 


(4.2.56) 


and 


^N(t*in)  £  K(Ctk)^fj  ^N(Wn-*-j)  ^  K(Ctk+j) 


I'o^Or^ 


;  j  =  0 


0  ;  else 


(4.2.57) 


For  large  M,  the  summation  of  equation  (4.2.55)  is  approximated  by  a  complex  Gaussian 

dksA 

random  variable  with  zero  mean  and  variance - Similarly,  the  term 

p(M+i)  ^ 


M 


-).[  f,)  ^N(^*^+m)  KlC^k+m) 


m  =  -M/2 


(4.2.58) 


is  approximated  by  a  complex  Gaussian  random  variable  with  zero  mean  and  variance 
p(M+l) 


Now  consider  the  tenn 


M/2 


L  Si  YN((On^n,)  ^  K^ttk  +m)  ' 


m=-M/2 


(4.2.59) 


and  note  that,  since  (ed))  and  |0(r)}  are  independent. 


68 


E  {^YN((0n)O*K(ak)  )  =  0  . 


(4.2.60) 


Furthermore, 


e|  YnCCOo)  ^  K(ttk)  J  ^  N  YN(Wn+j)  O  K(Ctk+j)  j 


IV  IV 

=  b(jjf  II  (ty  f  I  H«(Y)  H*, 


o  p  -Jt  -rc 

FN(0)n-(0o-  j)  F*N(Wn+j-C0p-  J)  FKCp-ttfc)  FKCp-ttk+j)  dy  d|3  .  (4.2.61) 


For  N  and  K  sufficiently  large. 


YN((0n)  0*K(ak)  1  f  YN(a)„,j)  0*K(ak.j) 


O^Oa^  2 

—^1  I  Ha(p(co„-(o„))  I  ;  j  =  0 

H  O 

0  ;  else 


(4.2.62) 


Therefd  e,  the  term  given  in  equation  (4.2.59)  is  approximated  by  a  zero  mean  complex 

CT^Oa^  2 

Gaussn  n  random  variable  with  variance - I  H„(p(ci)n-0)o))  I  . 

P(M+1)q 


Jur  estimate  for  the  variance  of  the  input  noise  process  is  given  by 


X  =  0 


(4.2.63) 


For  K  large,  we  may  approximate  b’  by 


(4.2.64) 


Our  estimate  for  the  system  transmission  function  is  now  aiven  bv 
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l^(O)n.aic)  =  S'M(0)n,ak)  , 


(4.2.65) 


or,  from  the  discussion  presented  above. 


a2+o^2 


Ho(p(a)n-(Oo))  +  ^Ao+)ll  +  ^2  +  il3  ;  tOn  =  COp  + 
P0+P1  +  P2+  P3  ;  else 


9^ 

p 


(4.2.66) 


where  po.  Pi.  P2.  and  P3  are  complex  zero-mean  Gaussian  random  variables  with 
variances 


'  p(M  +  l)(a2+0^2)  ' 

^ 

p(M+i)(g2+a^2)  • 


and 


03  2 


p(M>i)(g2+o^,2)  ^ 


Ho(p(a)n-a)o)) 


(4.2.67) 


(4.2.68) 


(4.2.69) 


(4.2.70) 


respectively. 

The  results  obtained  above  for  frequency-smoothed  estimates  of  K(a),a)  may  be 
extended  to  time-averajzed  estimates  for  the  system  transmission  function.  To  be  specific, 
assume  that  we  obtain  L  observations  of  the  input  and  output  signals  of  length  K  and  N, 
resfxvtively.  If  the  L  observations  of  the  output  process  are  obtained  in  such  a  manner  that 
the  input  samples  producing  one  set  of  output  samples  are  unrelated  to  the  input  samples 
prexlucing  any  other  set  of  output  samples  then  we  may  invoke  the  mixing  assumption  of 
Section  3.1  to  claim  that  the  L  output  observations  are  statistically  independent.  For  each 
observation,  form  the  discrete  F  ourier  transfonn  of  the  input  and  output  signal  in 
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accordance  with  equations  (4.2.6)  and  (4.2.11).  The  estimate  for  the  system  transmission 
function  is  obtained  from 

I^(o)„,ak)  =  SL((On,ak)  ,  (4.2.71) 

o2 

where  is  the  estimate  for  the  input  noise  variance  based  on  the  LK  available  samples  of 
the  input  process,  and 


SL(^n>Ctk) 


(4.2.72) 


In  equation  (4.2.72),  YN,i(0)n)  and  £K.i(«k)  are  the  1-th  di.screte  Fourier  transform  of  the 

input  and  output  signals.  Following  a  line  of  argument  identical  to  that  presented  for 
frequency-smoothed  estimates,  we  obtain 


I^(COn, 


CXk)  =  i 


2  Ho(p(a)n-COo))  +  Vo+V,-(-V2  +  V3  ;  (On  =  COo  -H 


Vq+V 1 +V2+ V3 


else 


oyc 

P 


(4.2.73) 


where  Vq.  Vj,  V2.  and  V3  are  complex  zero-mean  Gaussian  random  variables  with 
variances 


0,2  = 


'  pL(a2+aj,2)  ^ 


X  I  Ho(p(a)n-a)o)) 


2rf,2 


0^2  = 


0^0 


pL(02-t-0^2 


<720^2 

pL(a2-K0^2)  ’ 


and 


o 


4 


pL(a2-t-o^2) 


I  H„(p((0n-(0o))  r 


(4.2.74) 


(4.2.75) 


(4.2.76) 


(4.2.77) 
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respectively. 


Several  comments  can  be  made  regarding  the  method  for  estimating  K(a),a) 
presented  above.  First,  we  again  assumed  that  a  complex  input  signal  was  applied  to  the 
system  under  consideration.  For  systems  in  which  only  real  signals  may  be  applied  to  the 
system,  we  could  apply  the  real  and  imaginary  components  of  the  input  signal  to  the  system 
independently  and  linearly  combine  the  observed  outputs  to  estimate  the  system 
transmission  function,  as  was  described  in  Section  4.1.  However,  for  a  complex  white 
noise  process,  it  is  required  that  the  real  and  imaginary  components  be  independent  and 
identically  distributed  (see  Koopmans,  p.  263).  The  real  and  imaginary  components  for 
the  input  signal  could  therefore  be  obtained  from  two  K-lengih  observations  of  data 
obtained  from  a  single  white  noise  generator. 


Secondly,  we  assumed  that  the  input  noise  process  had  a  flat  spectrum,  i.e.,  we 
assumed  that  its  power  spectrum  was  constant  for  all  frequencies.  In  practice,  we  cannot 
guarantee  that  the  input  noise  process  has  a  flat  spectrum.  However,  we  could  estimate  the 
spectrum  of  the  input  noise  process,  Se(a),  from  the  given  observations,  and  estimate  our 
system  transmission  function  by 


A 


K(0)n,CXk) 


SMCtOnittk) 

Se(ak) 


(4.2.78) 
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4.3  CHAPTER  SUMMARY 


In  this  chapter  we  presented  two  techniques  for  empirically  characterizing  linear 
periodically  time  varying  systems.  The  properties  and  limitations  of  both  techniques  were 
discussed.  We  now  turn  our  attention  to  the  results  of  implementing  these  two  techniques 
to  the  characterization  of  several  LPTV  systems  of  interest. 
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5.  RESULTS  OF  IMPLEMENTATION  OF  METHODS  FOR  THE 
CHARACTERIZATION  OF  LINEAR 
TIME- VARYING  SYSTEMS 


We  now  briefly  present  the  results  of  implementation  of  the  methods  for 
characterizing  linear  periodically  time-varying  systems  presented  in  Chapter  4. 
Characterizations  of  a  set  of  LPTV  systems  were  conducted  using  both  deterministic  and 
nondeterministic  input  signals.  We  begin  by  considering  the  theoretical  system 
transmission  functions  for  the  LPTV  systems  of  interest. 

5.1  THEORETICAL  SYSTEM  TRANSMISSION  FUNCTIONS  FOR  THE 
LPTV  SYSTEMS  OF  INTEREST 

A  set  of  five  LPTV  systems  were  used  to  evaluate  the  methods  for  characterizing 
LPTV  systems  described  earlier.  The  first  two  of  these  five  systems  were  selected  because 
they  represented  basic  elements  that  are  commonly  found  in  more  complicated  LPTV 
systems.  Characterization  of  these  basic  system  elements  is  therefore  a  necessary  first  step 
towards  characterizing  more  complicated  LPTV  systems.  The  final  three  LPTV  systems 
evaluated  here  are  complicated  systems  based  on  the  interconnection  of  simpler  system 
elements. 


Figure  5.1.1  presents  the  block  diagram  of  the  first  LPTV  system  of  interest,  a 
simple  rate  1/2  decimator.  In  this  system,  the  input  signal,  x(T),  is  sampled  at  output 
sample  times,  t,  such  that  T  =  2t  to  produce  the  output  signal  y(l).  The  system  response 
function  ki(t.x),  is  therefore  given  by 


T  =  2t 
else 


(5.1.1) 


Using  equation  (2.6.8),  the  system  transmission  function  is  given  by 


Ki(co,a) 


-i(co-2a)t 

e 


-27tn-2a)  . 


(5.1.2) 


oo 


OO 


This  function  is  indicated  schematically  in  Figure  5.1.2.  For  this  figure,  the  squared 
magnitude  of  Ki(a),a)  was  evaluated  at  the  set  of  discrete  points  in  the  ((0,a)  bifrequency 
plane  defined  by 


27tn 

C0n  =  -J^  ,  (5.1.3) 

and 

27tk 

ak  =  -|^  (5.1.4) 

where  N  =  128  and  K  =  64.  Note  that  input  and  output  frequencies  have  been 
nonnalized  by  k  in  Figure  5.1.2. 

Figure  5.1.3  pre.sents  the  block  diagram  for  the  second  LPTV  system  of  interest,  a 
simple  quadrature  demodulator.  In  this  system,  the  input  signal,  x(t),  is  demodulated  by  a 
complex  sinusoid  of  frequency  cUo  to  obtain  the  output  signal  y(T).  The  system  response 
function,  k2(t,x),  is  therefore  given  by 

I  -nco-ajo)i 

Mt.i:)=  •  ^5.1.5) 

I  0  ;  else 


The  system  transmission  function  for  this  sy.stem  is  given  by 

oo 

K2(0),a)  =  c 

I  =  -  oo 


(5.1.6) 


or 


K;(CO.(X) 


I 


0((o-2;in+to„-a) 


n  =  CO 


(5.1.7) 


I  his  svstem  transmission  function  is  indicated  schematically  in  Figure  5.1.4.  The  squared 
magnitude  of  K2(0).a)  was  again  evaluated  at  a  discrete  set  of  points  in  the  (oo.a) 


bifrequency  plane  according  to  equations  (5.1.3)  and  (5.1.4)  with  N  =  K  =  64.  Input 
and  output  frequencies  have  been  normalized  by  k  in  Figure  5.1.4. 

Figure  5.1.5  presents  the  block  diagram  for  the  third  LPTV  system  of  interest,  a 
rate  3/2  interpolation  system.  In  this  system,  the  input  signal  is  up-sampled  (i.e.,  zero- 
padded)  by  a  factor  of  three,  passed  through  a  LTI  low  pass  filter,  and  decimated  by  a 
factor  of  two.  The  system  response  function  for  the  rate  3  up-sampling  is  given  by 

k(t.T)={  •  (5.1.8) 

I  0  ;  else 

The  system  transmission  function  for  the  rate  3  up-sampling  is  therefore  given  by 


K^p,a) 


5(P  2Tin-a/3)  . 

n  =  —  po 


(5.1.9) 


The  overall  system  transmission  function  for  the  rate  .3/2  interpolation  system,  Ko(C0,a),  is 
given  by 


Ko(co,a) 


K,(o),p)  H(())  K(p.a)dp  , 


(5.1.10) 


or  simply 


Ko(w.a) 


yo  oc 


2  H(a/3+27tn)  6(to-27i(m-2n)-2a/3) 


m  =  -  «>  n  =  -  o<’ 


(5.1.11) 


Figure  5. 1.(1  schematically  indicates  this  system  transmission  function.  Figure  5.1.7 
presents  a  three-dimensional  view  of  K<,((o,a).  For  the.se  plots,  the  .souared  magnitude  of 
the  system  transmission  function  was  evaluated  at  a  .set  of  discrete  points  in  the  (w,a) 
bifrequency  plane  according  to  equations  (5.1.3)  and  (5.1.4)  with  N=96  and  K=64. 

Figure  5.1.8  presents  a  block  diagram  for  the  fourth  LPTV  system  of  interest,  a 
rate  1/2  quadrature  sampling  system.  In  this  system,  the  input  signal,  x(T),  is  demodulated 
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by  cos  x)  and  -i  sin  (^x),  and  sampled  in  quadrature  to  obtain  the  in-phase  and 
quadrature  components  of  a  digital  quadrature  representation  for  x(x)  (see  Baugh,  1988). 
The  system  response  function  for  this  system  is  given  by 


k(t,x)  =  i 


COS(2  X  ) 


sin(^  X ) 
0 


X  =  2 1 

X  =  2t+ 1 
e !  se 


(5.1.12) 


The  system  transmission  function  for  the  rate  1/2  quadrature  sampling  system  is  therefore 
given  by 


K(to,a)  =  (I 


■i  e‘“) 


8{a)-27tn-2a-ii)  +  6(ca-27tn-2a-t-7i) 


(5.1.13) 


n  =  - 


This  system  transmission  function  is  indicated  schematically  in  Figure  5.1.9,  with  the 
input  and  output  frequencies  again  normalized  by  7i.  Figure  5.1.10  presents  a  three- 
dimensional  view  of  this  system  transmi.ssion  function.  Again,  the  squared  magnitude  of 
the  system  transmission  function  was  evaluated  at  a  set  of  discrete  points  in  the  (co,a) 
bifrequency  plane  according  to  equations  (5.1.3)  and  (5.1.4)  with  N=128  and  K=64. 
Figure  5.1.11  presents  the  result  of  plotting  the  largest  value  of  K(co,a)  over  all  O)  for  each 
a.  In  a  sense,  this  plot  gives  an  indication  of  the  response  of  the  LPTV  system  to  each 
input  frequency.  All  frequencies  have  been  n^miabzed  by  it.  Consider  that  on  the  lines 

0)  =  2a-t-7t  (5.1.14) 


and 


to  =  2a-n 


(5.1.15) 


the  .square-magnitude  of  K(a),a)  is  given  by 

I  K(to,a)  i’ =  2(  1 +sin  a)  .  (5.1.16) 

This  function  is  clearly  indicated  in  Figure  5.1.1 1. 
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Figure  5.1.12  presents  the  block  diagram  for  the  fifth  LPTV  system  of  interest,  a 
rate  1/2  quadrature  demodulator.  Note  that  this  system  is  a  series  connection  of  the  simple 
quadrature  demodulator  of  Figure  5.1.3  with  a  gain  of  2,  a  linear  time-invariant  low  pass 
filter,  and  the  rate  1/2  decimator  of  Figure  5.1.1.  Using  the  results  of  Section  2.4  for  the 
analysis  of  LTV  network  structures,  we  may  express  the  system  transmission  function  for 
‘.his  LPfV  system,  K3(a),a),  as 


K3(co.a)  =  ^  J  2Ki((o,l3)  H(p)  KiCP.a)  d|3  , 


(5.1.17) 


where  H(p)  is  the  frequency  response  of  the  LTI  system.  Substituting  equations  (5.1.2) 
and  (5.1.7)  into  (5.1.17),  we  obtain 


K3(C0 


,«)=  I  I 


2  H(a-c0o+2nn)  5({o-2Ti(m-2n)+2o)o-2a)  .  (5.1.18) 


m  =  -oon  =  -oo 


This  system  transmission  function  is  indicated  schematically  in  Figure  5.1.13,  with  the 
input  and  output  frequencies  again  normalized  by  7t.  Figure  5.1.14  presents  a  three- 
dimensional  view  of  this  system  transmission  function,  indicating  the  modulation  of  the 
impulse  sheet  by  H(P),  the  LTI  frequency  response  function.  For  these  plots,  the  squared 
magnitude  of  the  system  transmission  function  was  evaluated  at  a  set  of  discrete  points  in 
the  (CO, a)  bifrequency  plane  according  to  equations  (5.1.3)  and  (5.1.4)  with  N=128  and 
K=64.  Figure  5.1.15  presents  the  result  of  plotting  the  largest  value  of  K3(to,a)  over  all 
(0  for  each  a. 
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Output  Frequency 
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AS-9 1-371 

FIGURE  5  13 

BLOCK  DIAGRAM  FOR  THE  SIMPLE  QUADRATURE  DEMODULATOR 


Input  Frequency 
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FIGURE  5.L4 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR 
A  SIMPLE  QUADRATURE  DEMODULATOR 
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FIGURE  5  15 
BLOCK  DIAGRAM  FOR  THE  RATE  3/2  INTERPOLATOR 
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FIGURE  5.1.6 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE 
RATE  3/2  INTERPOLATION  SYSTEM 
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FIGURE  5.1.9 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE 
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FIGURE  5.1.10 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE 
RATE  1/2  QUADRATURE  SAMPLINi;  SYSTEM 
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FIGURE  5.1.11 

THEORETICAL  SYSTEM  RESPONSE  FOR  THE  RATE  1/2 
QUADRATURE  SAMPLING  SYSTEM 
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FIGURE  5.1.17 

BLOCK  DIAGRAM  FOR  THE  RATE  1/2  QUADRATURE  DEMODULATOR 
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FIGURE  5.1.13 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE 
RATE  1/2  QUADRATURE  DEMODULATION  SYSTEM 
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FIGURE  5.1.14 

THEORETICAL  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE 
RATE  1/2  QUADRATURE  DEMODULATION  SYSTEM 
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FIGURE  5.1.15 

THEORETICAL  SYSTEM  RESPONSE  FOR  THE  RATE  1/2 
QUADRATURE  DEMODULATION  SYSTEM 
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5.2  CHARACTERIZATION  via  DETERMINISTIC  SIGNALS 


The  method  of  characterizing  LPTV  systems  using  deterministic  signals  described 
in  Section  4. 1  was  applied  to  the  five  systems  of  interest.  A  complex  sinusoid  of  a 
specific  frequency  was  generated  and  applied  to  the  LPTV  system  of  interest.  A 
512-sample  record  of  the  system  output  was  obtained,  and  a  minimum  3-term  Blackman- 
Harris  windowing  function  was  then  applied  to  this  output  record  (see  Harris,  1978).  The 
512-point  discrete  Fourier  transform  of  this  windowed  output  data  record  was  computed. 
The  resulting  DFT  was  smoothed  using  M  =  8  in  order  to  reduce  the  number  of  points  in 
the  resulting  estimate  for  K(ct),a),  and  so  simplify  the  plotting  of  the  data.  This  64-point 
smoothed  DFT  was  used  as  the  estimate  for  the  system  transmis  ion  function  for  the 
specific  input  frequency.  The  procedure  was  then  repeated  for  each  of  K  discrete  input 
frequencies,  Uk  =  Y’  where  K  was  chosen  so  that  N/K  =  p,  the  rate  of  the  LPTV  system 
under  evaluation. 

Figures  5.2.1  and  5.2.2  present  the  estimated  system  transmission  functions  for 
the  simple  rate  1/2  decimator  and  the  simple  quadrature  demodulator  obtained  using  this 
method.  Note  the  close  agreement  between  these  estimates  and  the  theoretical  system 
transmission  functions  of  Figures  5.1.2  and  5.1.4.  Figure  5.2.3  presents  the  estimated 
system  transmission  function  for  the  rate  3/2  interpolator  as  a  logarithmic  contour  plot. 
Figure  5.2.4  presents  a  three-dimensional  view  of  this  estimated  system  transmission 
function.  Figure  5.2.5  presents  the  estimated  system  transmission  function  for  the 
rate  1/2  quadrature  sampling  system.  Figure  5.2.6  presents  a  three-dimensional  view  of 
the  estimated  transmission  function.  Figure  5.2.7  presents  the  result  of  plotting  the  largest 
value  of  the  estimate  for  Klco.a)  over  all  (O  for  each  a.  Note  the  close  agreement  between 
the  estimated  system  transmission  function  and  the  theoretical  value  plotted  in 
Figures  5. 1 .9,  5.1.10  and  5.1.11.  Figure  5.2.8  presents  the  estimated  system 
transmission  function  for  the  rate  1/2  quadrature  demodulator.  In  this  contour  plot,  the 
magnitudes  are  presented  in  logarithmic  form.  Figure  5.2.9  presents  a  three-dimensional 
view  of  the  estimated  system  transmission  function  for  the  rate  1/2  quadrature 
demodulator.  Figure  5.2. 10  presents  the  result  of  plotting  the  largest  value  of  K(to,a)  over 
all  0)  for  each  a.  Again  there  is  clo.se  agreement  between  the  estimated  system  transmission 
function  and  the  theoretical  values  presented  in  Figures  5.1.13,  5.1.14,  and  5.1.15. 

A  simulation  of  the  effects  of  quadrature  phase  error  between  the  real  and  imaginary 
components  of  the  input  reference  sinusoid  was  also  conducted.  Specifically,  the  real  and 
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imaginary  parts  of  the  input  reference  sinusoid  were  generated  and  applied  separately  to  the 
input  of  the  rate  1/2  quadrature  demodulator.  A  quadrature  phase  error  of  1°  was 
introduced  between  the  real  and  imaginary  components  of  the  input  reference  sinusoid. 
From  the  discussion  of  Section  4.1,  we  would  expect  to  see  the  term  |  K(a),-a)  appear  in 

our  estimate,  where  |  is  roughly  -41  dB.  Figure  5.2.1 1  presents  a  contour  plot  of  the 

logarithm  of  the  estimated  system  transmission  function  obtained  using  the  given  reference 
sinusoids.  Figure  5.2.12  presents  the  result  of  plotting  the  largest  value  of  K(a),a)  over 
all  0)  for  each  a.  From  these  two  figures  we  note  that  quadrature  phase  error  in  the  input 
reference  sinusoids  does  indeed  lead  to  the  term  |  K(o),-a)  in  our  estimated  system 
transmission  function. 
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FIGURE  5.2.5 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  SAMPLING  SYSTEM 


AS-91-389 


FIGURE  5.2.6 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  SAMPLINC;  SYSTEM 
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FIGURE  5.2.7 
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FIGURE  5.2.8 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 

QUADRATURE  DEMODULATOR 


System  Response  (dB) 


AS-9 1-392 

FIGURE  5.2.9 
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FIGURE  5.2.10 
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FIGURE  5.2.11 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  DEMODULATION  SYSTEM  WITH  1  DEGREE  OF 
QUADRATURE  PHASE  ERROR 
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5.3  CHARACTERIZATION  via  NONDETERMINISTIC  SIGNALS 


The  method  of  characterizing  LPTV  systems  using  nondeterministic  input  signals 
was  used  to  estimate  the  system  transmission  functions  for  the  five  LF*TV  systems  of 
interest.  In  each  case,  a  Gaussian  input  sequence  was  applied  to  the  system  of  interest,  and 
an  estimate  of  the  system  transmission  function  was  obtained  by  the  appropriate  time  and 
frequency  smoothing  of  the  discrete  spectral  components.  Normalization  of  the  estimated 
system  transmission  function  by  the  estimated  input  noise  power  spectrum  was  used,  as 
described  in  Section  4.2. 

Figure  5.3.1  presents  a  contour  plot  of  the  logarithm  of  the  estimated  system 
transmission  function  for  the  simple  rate  1/2  decimator.  For  this  simulation,  the  input  DFT 
length,  K,  was  512,  the  output  DFT  length,  N,  was  1024,  the  frequency  smoothing  length, 
M,  was  8,  and  the  number  of  averages,  L,  was  512.  Figure  5.3.2  presents  the  contour 
plot  of  the  logarithm  of  the  estimated  system  transmission  function  for  the  simple 
quadrature  demodulator.  For  this  simulation,  K=N=512,  M=8,  and  L=512. 

Figure  5.3.3  presents  a  contour  plot  of  the  logarithm  of  the  estimated  system 
transmission  function  for  the  rate  3/2  interpolator.  Figure  5.3.4  presents  a  three- 
dimensional  view  of  this  estimate.  For  this  simulation,  K=512,  N=768,  M=8,  and  L=512. 

Figure  5.3.5  presents  a  contour  plot  of  the  logarithm  of  the  estimated  system 
transmission  function  for  the  rate  1/2  quadrature  sampling  system.  Figure  5.3.6  presents 
a  three-dimensional  view  of  this  estimate.  Figure  5.3.7  presents  the  result  of  plotting  the 
largest  value  of  K(a),a)  over  all  (o  for  each  a.  For  this  simulation,  K=512,  N=1024, 
M=8.  and  L=512. 


Figure  5.3.8  presents  a  contour  plot  of  the  logarithm  of  the  estimated  system 
transmission  function  for  the  rate  1/2  quadrature  demodulator.  Figure  5.3.9  presents  a 
three-dimensional  view  of  this  estimate.  Figure  5.3.10  presents  the  result  of  plotting  the 
largest  value  of  K((o,a)  over  all  w  for  each  a.  For  this  simulation,  K=512,  N=1024, 
M=8,  and  L=512.  Note  that  we  obtain  rea.sonably  good  agreement  with  the  theoretical 
response  presented  in  Section  5.1. 

A  simulation  of  the  effects  of  using  separate  real  and  imaginary  input  signal 
components  was  also  conducted.  For  this  simulation,  a  single  Gaussian  noi.se  .sequence 
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was  generated  and  applied  to  the  input  of  the  rate  1/2  quadrature  demodulation  system. 
Alternating  records  of  the  observed  input  and  output  sequences  were  used  to  produce  the 
real  and  imaginary  parts  of  the  estimated  system  transmission  function.  As  in  the  previous 
simulation  for  the  rate  1/2  quadrature  demodulation  system,  K=512,  N=1024,  M=8,  and 
L=512.  Figure  5.3.1 1  presents  the  contour  plot  of  the  logarithm  of  the  estimated  system 
transmission  function  obtained  from  this  simulation.  Figure  5.3.12  presents  a  three- 
dimensional  view  of  the  estimate. 


98 


Output  Frequency 


Input  Frequency 


AS-91-396 

FIGURE  5.3.1 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  SIMPLE 

RATE  1/2  DECIMATOR 


Input  Frequency 

AS-9 1-397 

FIGURE  5.3.2 

ESTIMATED  SYS  i  EM  TRANSMISSION  FUNCTION  FOR  THE  SIMPLE 

QUADRATURE  DEMODULATOR 


99 


Output  Frequency 


Input  Frequency 


AS-91-398 
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FIGURE  5.3.5 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  SAMPLING  SYSTEM 
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FIGURE  5.3.6 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  SAMPLING  SYSTEM 
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FIGURE  5.3.9 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
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FIGURE  5.3.10 

ESTIMATED  SYSTEM  RESPONSE  FOR  THE  RATE  1/2  QUADRATURE 

DEMODULATOR 


103 


Output  Frequency 


Input  Frequency 


AS-9 1-406 

FIGURE  5.3.11 

ESTIMATED  SYSTEM  TRANSMISSION  FUNCTION  FOR  THE  RATE  1/2 
QUADRATURE  DEMODULATION  SYSTEM  USING  REAL  INPUT 
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FIGURE  5.3.12 
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5.4  CHAPTER  SUMMARY 


In  this  chapter,  the  techniques  for  characterizing  LPTV  systems  described  in 
Chapter  4  were  applied  to  several  systems  of  interest.  The  consequences  of  quadrature 
phase  error  and  independent  application  of  the  real  and  imaginary  components  of  the  input 
signal  were  also  demonstrated. 
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6.  CONCLUSIONS 


In  this  paper  we  have  presented  two  techniques  for  the  empirical  characterization  of 
discrete-time  linear  periodically  time-varying  systems.  These  techniques  were  developed 
amid  a  framework  of  concepts  related  to  the  analysis  of  discrete-time  LPTV  systems  and  to 
the  signals  produced  by  such  systems.  The  results  of  applying  these  techniques  to  several 
LPTV  systems  of  interest  were  also  presented. 

From  the  tests  presented  in  Chapter  5,  we  see  that  characterization  of  LPTV 
systems  using  deterministic  sinusoidal  reference  inputs  provides  an  accurate  estimate  of  the 
system  transmission  function  only  when  the  quadrature  phase  error  between  the  real  and 
imaginary  components  of  the  reference  input  is  maintained  within  some  tight  limit  over  the 
analysis  bandwidth  of  interest.  In  those  situations  where  we  may  apply  only  real  input 
signals,  so  that  the  real  and  imaginary  components  must  be  applied  separately,  this 
requirement  may  place  a  severe  limitation  on  the  design  of  the  reference  signal  generator. 
In  contrast,  the  use  of  non-deterministic  input  signals  for  characterization  of  LPTV  systems 
may  proceed  by  applying  independent  real  and  imaginary  input  signal  components  to  the 
system  of  interest.  The  design  of  the  reference  signal  generator  is  therefore  greatly 
simplified. 

In  conclusion,  the  characterization  of  discrete-time  LPTV  systems  using 
deterministic  signals  is  appropriate  in  those  situations  where  a  fiiie  estimate  of  the  system 
transmission  function  is  required  at  a  particular  frequency,  or  across  some  relatively  narrow 
band  of  frequencies.  The  use  of  non-deterministic  signals  for  LPTV  system 
characterization  is  appropriate  in  those  situations  where  an  estimate  of  the  system 
transmission  function  is  required  over  a  broad  range  of  frequencies,  or  where  the  direct 
application  of  a  reference  sinusoid  would  be  impractical. 
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25  Library 


Commander 

Naval  Space  and  Warfare  Systems  Command 
Department  of  the  Navy 
Washington,  D.C.  20.363-5101 

26  Attn:  Library 

Commanding  Officer 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
Stennis  Space  Center,  MS  39529-5001 

27  Attn:  Library 

Commanding  Officer 

Naval  Oceanographic  Office 

Stennis  Space  Center,  MS  39522-5(X)l 

28  Attn:  Library 

Director 

Naval  Research  Laboratory 
Washington,  D.C.  20375-5(KX) 

29  Attn:  Library 

Officer  in  Charge 
Naval  Surface  Warfare  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 

30  Attn:  Library 
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31 

32 

33 

34 

35 

36 

37-39 

40 

41 


Commanding  Officer 
Naval  Coastal  Systems  Center 
Panama  City,  FL  32401 
Attn:  Library 

Officer  in  Charge 
Naval  Ocean  Systems  Center 
Hawaii  Laboratory 
Kailua,  HI  96734-0977 
Attn:  Library 

Office  of  the  Chief  of  Naval  Research 
Office  of  Naval  Technology 
Department  of  the  Navy 
Arlington,  VA  22217-5000 
Attn:  Library 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington.  VA  22209-2309 
Atui:  Library 

National  Aeronautics  and  Space  AdminisU'aiion 
Langley  Research  Center 
Hampton,  VA  23665 
Attn:  Library 

Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
Washington,  D.C.  20315-0030 
Attn:  Library 

Commanding  Officer  and  Director 
Defense  Technical  Information  Center 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria.  VA  22314 

Tracor  Applied  Sciences,  Inc. 

AARD/Sca  Operations 
35  Thomas  Griffin  Rd. 

New  London,  CT  06320 
Attn:  D.  Turbes 

Planning  Systems,  Inc. 

792  S.  Westpark  Dr. 

McLean,  VA  22102 
Attn:  H.  Wiu'ner 
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42 

Lockheed  Missile  and  Space  Corporation 
0/6580,  B/582N 

1 1 1 1  L(x:kheed  Way 

Sunnyvale,  CA  94088-3504 

Attn:  W.  Eugene  Brown 

43 

National  Academy  of  Sciences 

National  Research  Council 

2101  Constitution  Avenue 

Washington,  D.C.  20418 

Attn:  Library 

44 

Applied  Research  Laboratory 

The  Pennsylvania  State  University 

P.O.  Box  30 

State  College,  PA  16801 

Attn:  Library 

45 

Applied  Physics  Laboratory 

The  Johns  Hopkins  University 

Johns  Hopkins  Road 

Laurel,  MD  20810 

Attn:  Library 

46 

Applied  Physics  Laboratory 

The  University  of  Washington 

1013  Northeast  40th  Street 

Seattle,  WA  98015 

Attn:  Library 

47 

Marine  Physical  Laboratory 

The  Scripps  Institution  of  Oceanography 
The  University  of  California  at  San  Diego 
San  Diego,  CA  92152 

Attn:  Library 

48 

Woods  Hole  Oceanographic  Institution 
86-95  Water  Street 

Woods  Hole,  MA  02543 

Attn:  Library 

49 

Signal  Physics  Group,  ARL:UT 

50 

Martin  L.  Barlett,  ARL:UT 

51 

Garland  R.  Barnard,  ARL.UT 

52 

Kevin  W.  Baugh,  ARL:UT 

53 

M.  Gette  Beauvais,  ARL.UT 
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54 

Glen  E.  Ellis,  ARL;UT 

55 

Joseph  F.  England,  ARL:UT 

56 

Douglas  J.  Fox,  ARLrUT 

57 

Marshall  E.  Frazer,  ARL:UT 

58 

Keith  R.  Hardwicke,  ARL;UT 

59 

Melvin  J.  Hinich,  ARL:UT 

60 

James  L.  Lamb,  ARL.UT 

61 

Fredrick  W.  Machell,  ARL:UT 

62 

Chetna  S.  Mathur,  ARL.UT 

63 

Sara  S.  Matzner,  ARL:UT 

64 

Thomas  G.  Muir,  ARL;UT 

65 

Clark  S.  Penrod,  ARL.UT 

66 

F.  xMichacl  Pestorius,  ARL:UT 

67 

Edward  J.  Powers,  ARLrUT 

68 

Jack  H.  Sheehan,  ARLrUT 

69 

Ronald  O.  Steannan,  ARLrUT 

70 

Robert  T.  Trochta,  ARLrUT 

71 

J.  Kenneth  Vaughan,  ARLrUT 

72 

Joseph  F.  Willinan,  ARLrUT 

73 

Gary  R.  Wil.son,  ARLrUT 

74 

Library,  ARLrUT 

75  -  85 

Reserve,  ARLrUT 
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